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Abstract. This talk is an attempt to clarify for experimentalists what do we (the¬ 
orists) mean when we are saying “Pomeron”. I hope, that in this talk they will find 
answers to such questions as what is Pomeron, what we have learned about Pomeron 
both experimentally and theoretically, what is the correct strategy to study Pomeron 
experimentally and etc. I also hope that this talk could be used as a guide in the zoo 
of Pomerons: “soft” Pomeron, “hard” Pomeron, the BFKL Pomeron, the Donnachie - 
Landschoff Pomeron and so on. The large number of different Pomerons just reflects 
our poor understanding of the high energy asymptotic in our microscopic theory - 
QCD. The motto of my talk, which gives you my opinion on the subject in short, is: 
Pomeron is still unknown but needed for 25 years to describe experimen¬ 
tal data . This is a status report of our ideas, hopes, theoretical approaches and 
phenomenological successes that have been developed to achieve a theoretical under¬ 
standing of high energy behaviour of scattering amplitude. 


I HIGH ENERGY GLOSSARY 

In this section I am going to define terminology which we use discussing high 
energy processes. Most of this terminology has deep roots in Reggeon approach 
to high energy scattering and it will be more understandable after my answer to 
the question:“what is the “soft” Pomeron”. Although I cannot explain right now 
everything, I think, it will be very instructive to have a high energy glossary in 
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FIGURE 1. The high energy interaction in the parton model. 


front of your eyes since it gives you some feeling about problems and ideas which 
are typical for high energy asymptotic. 

To make easier your understanding I give here the picture of a high energy 
interaction in the parton model ( see Fig.l ). 

In the parton approach the fast hadron decays into point-like particles ( partons ) 
a long before ( typical time i oc ^ ) the interaction with the target. However, 
during time t all these partons are in a coherent state which can be described by 
means of a wave function. The interaction of the slowest parton ( “wee” parton ) 
with the target destroys completely the coherence of the partonic wave function. 
The total cross section of such an interaction is equal to 

cr tot = N x (To (1) 


where 

• N = flux of “wee” partons ; 

• cr 0 = the cross section of the interaction of one “wee” parton with the target . 

One can see directly from Fig.l that the flux of “wee” partons is rather large 
since each parton in the parton cascade can decay into its own chain of partons. 
If one knows the typical multiplicity < n > of a single chain ( say, dotted one in 
Fig.l ) he can evaluate the number of “wee” partons. Indeed, 

N = e <n> (2) 

where < n > oc ^ and Ay is the average distance in rapidity for partons in one 
chain. Reasonable estimate for Ay is Ay ~ A-. Therefore, we have finally 

< n > — lj 0 Ins ; u 0 = Const a s , (3) 

It is obvious that Eq. (3) leads to 


IV oc s“V 












Finally, Eq. (1) can be rewritten in the form: 

§ 

a t ot = P(projectile) x ao(target) x ( — , (4) 

so 

where P is a coefficient which depends only on the quantum numbers and variables 
related to the projectile and which is a physics meaning of probability to find 
the parton cascade of Fig.l in the fast hadron (projectile). However, our picture 
should be relativistic invariant and independent of the frame in which we discuss 
the process of scattering. Therefore, P in Eq. (4) have to be proportional to 
P oc a °(P r °P ctlle ) an( ] ( 4 ) i s be reduced to 

1 / s \ UJ ° 

cr to t = cr 0 (projectile) x o 0 (target) x — x — , (5) 

So V So J 

where So is the energy scale which says what energy we can consider as high one 
for the process of interaction. High energy asymptotic means s > «o and namely 
the value of s 0 is the less known scale in high energy physics. 

ft should be stressed that the value of u 0 does not depend on any variable related 
both to projectile and/or target. It only depends on density of partons in the 
partons cascade or, in other words, on the parton emission in our microscopic 
theory. 

1. Reggeon: 

The formal definition of Reggeon is the pole in the partial wave in /-ch ann el of 
the scattering process. For example, for n + + p clastic amplitude the Reggeon 
is a pole in the partial wave of the reaction: 7 r + + 7 r~ —> p + p, namely, the 
amplitude of this process can be written in the form: 

OO 

m = E fi(t)(2l + l)P,(z) , (6) 

1=0 

where z = cosO and 6 is the scattering angle from initial pion to final proton 
(antiproton). Reggeon is the hypothesis that fi(t) has a pole of the form 

«» - 

where function a#(t) is the Reggeon trajectory which experimentally has a 
form: 


a R (t) = a R ( 0 ) + a' R (0)t , 


( 8 ) 


where 0:^(0 is the intercept of the Reggeon and «/,(()) is its slope. It turns out 
that such hypothesis gives the following asymptotic at high energy ( see Refs. 
[ 1 ] [ 2 ] ) 



s «rM ( — g\ a R(f) 

A P (s, t) = gi(mi, Mi, t) x g 2 (m 2 ,M 2 ,t) x - : - — - . (9) 

svmi a r{i) 

This equation has two important properties: (i) at t — m 2 R , where rn R is the 
mass of resonance with spin j ( j = a R {t = m R ) ) it describe the exchange of 
the resonance, namely, 

S’ j 

A R (s,t -*• m|) -> gig 2 — -- 

m R — t 

and (ii) in the scattering kinematic region ( t < 0 ) A(s,t ) oc which 

gives the asymptotic for the total cross section 

<7 tot oc < Coust in 2 s(Froissart bound) 

for or( 0) < 0. Summarizing what I have discussed I would like to emphasize 
that the Reggeon hypothesis is the only one which has solved the contradiction 
between experimental observation of resonances with spin larger than 1 and a 
violation of the Froissart theorem which, at first sight, the exchange of such 
resonances leads to. 

2. “Soft” Pomeron: 

Actually, “soft” Pomeron is the only one on the market which can be called 
Pomeron since the definition of Pomeron is 

Pomeron is a Reggeon with the intercept close to 1, or in other 
words, o:p(0) — 1 = A < 1. 

1 want to point out three extremely important experimental facts, which were 
the reason of the Pomeron hypothesis: 

1. There are no resonances on the Pomeron trajectory; 

2. The measured total cross sections are approximately constant at high en¬ 
ergy; 

3. The intercepts of all Reggeons are smaller than 1 and, therefore, an ex¬ 
change of the Reggeon leads to cross sections which fall down at high energies 
in the contradiction with the experimental data. 

It should be stressed that there is no deep theory reason for the Pomeron 
hypothesis and it looks only as a simplest attempt to comply with the exper¬ 
imental data on total cross sections. The common belief is that the “soft” 
Pomeron exchange gives the correct high energy asymptotic for the processes 
which occur at long distances ( “soft” processes ). However, we will discuss 
below what is the theory status of our understanding. 

3. Pomeron structure: 

We say “Pomeron structure” when we would like to understand what inelastic 
processes are responsible for the Pomeron exchange in our microscopic theory. 



In some sense, the simplest Pomeron structure is shown in Fig.l in the parton 
model. Unfortunately, we have not reached a much deeper insight than it is 
given in the picture of Fig.l. 

4. Donnachie - Landshoff Pomeron: 

The phenomenology based on the Pomeron hypothesis turned out to be very 
successful. It survived at least two decades despite of a lack of theoretical 
arguments for Pomeron exchange. Donnachie and Landshoff [3] gave an el¬ 
egant and economic description almost all existing experimental data [3] [4] 
assuming the exchange of the Pomeron with the following parameters of its 
trajectory: 

1. A = ap(0) - 1 « 0.08 ; 

2. a' P (0) « 0.25 GeV~ 2 . 

The more detail properties of the D-L Pomeron we will discuss below. At 
the moment, we can use these parameters when we are going to make some 
estimates of the Pomeron contribution. 

5. “Hard” Pomeron: 

“Hard” Pomeron is a substitute for the following sentence: the asymptotic for 
the cross section at high energy for the “hard” processes which occur at small 
distances (r_i_) of the order of 1 /Q where Q is the largest transverse momentum 
scale in the process. The “hard” Pomeron is not universal and depends on 
the process. The main brick to calculate the “hard” Pomeron contribution is 
the solution to the evolution equation ( DGLAP evolution [5] ) in the region 
of high energy. The main properties of the “hard” Pomeron which allow us to 
differentiate it from the “soft” one are: 

1 . the energy behaviour is manifest itself through the variable xb = — ; 

2. the high energy asymptotic can be parameterized as a to t oc but the 

power uq(Q 2 ) crucially depends on Q 2 . Recall that it is not the case for the 
“soft” Pomeron or any Reggeon; 

3. although the energy and Q 2 behaviour of the “hard” Pomeron is quite 
different from the “soft” one the space - time picture for it is very similar to 
that one which is shown in Fig. 1. 

The procedure of calculation of the “hard” Pomeron is based on the perturba¬ 
tion theory. Indeed, for “hard” processes we have a natural scale of hardness 
that we call Q 2 3> m 2 where m 2 is the typical mass scale for “soft” inter¬ 
actions. The QCD coupling constant as(Q 2 ) can be considered as a small 
parameter ( a s (Q 2 ) <C 1 ) with respect to which we develop the perturba¬ 
tion theory. Generally speaking any physics observable ( let say the gluon 
structure function xbG(xb,Q 2 ) can be written as a perturbative series in the 
form: 



x b G(x b ,Q 2 ) = lim | Q2>>m 2 a7 s(Q 2 ) x M n (x B ,Q' 2 ) . (10) 

n =0 

For “hard” processes we change the order of summation and limit 1 and get 

OO 

x b G(x b ,Q 2 ) = a s(Q 2 ) x \q 2 > m 2 M n (x B , Q 2 ) ■ (11) 

n =0 


The analysis of the Feyman diagrams shows that for Q 2 S> m 2 M n (x,Q 2 ) 
has a form 


li m | 0 2» m 2M„(a; B , Q 2 ) => (12) 

Pn(x)L" + p£\x) L- 1 + ... + f>„” +1 (x) , 

where L is a large log in the problem, namely, L = In (Q 2 /m 2 ). 

Taking only the leading term with respect to L ( L n ) in Eq. (12), we arrive to 
so called leading log approximation (LLA) of perturbative QCD: 

OO 

x B G LLA (x Bl Q 2 ) = £ a n s {Q 2 ) x P n {x)L n ,. (13) 

n=0 


The practical way how we perform summation in Eq. (13) is the DGLAP 
evolution equation [5] and LLA is the same as the solution of the DGLAP 
evolution equation in the leading order ( LO DGLAP ). Taking a next term in 
L in Eq. (12) of the order of L n_1 in addition to the leading one we can arrive 
to the gluon structure function in the next to leading order (NLO DGLAP ). 
As I have mentioned, the DGLAP evolution equation gives a regular procedure 
for calculation . 

I think, that it is instructive to show here the “hard” Pomeron contribution to 
the gluon structure function in the region of large Q 2 Qq and small x B ( see 
Ref. [6] for example). It turns out that in this case P n {x) has a leading term of 
the order ln n (l/a;), reducing the problem to summation of perturbative series 
with respect to parameter ag ln(l/x) ln(Q 2 /m 2 ): 

x b G dla (x b ,Q 2 ) = ^2,n — 0°°c n (ag ln(l/x) ln(Q 2 /m 2 ) ) n . (14) 


The summation over n which we can also perform using the DGLAP evolution 
equation in so called double log approximation ( DLA) gives 


x b G dla (x b ,Q 2 ) 
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e 


2 N c 
b 0 


\n(l/xg) 111 


a s(Qp) 


(15) 


It should be stressed that this change is not so obvious but fortunately for this particular case 
we can use the powerful method of the renormalization group approach to justify it a posteriori. 



where N c is the number of colours and b 0 = (11 — |7Vy)/2 with Nf flavours. 
The running QCD coupling is equal to as(Q 2 ) = hQ . 

One can see from Eq. (15) that the real asymptotic or the “hard” Pomeron is 
quite different from the Reggeon - like behaviour and could be approximated 
by power - like function in very limited kinematic region of xi , and Q 2 . 1 believe, 
that this fact should be known to any experimentalist even if he is fitting the 
data using the Regge - like parameterization xbG(xb,Q 2 ) = const■ 

6. The BFKL Pomeron: 

The BFKL Pomeron is what we have in perturbative QCD instead of “soft” 
Pomeron. The method of derivation for the BFKL Pomeron is designed spe¬ 
cially for the processes with one but rather large scale ( Q 2 ~ Ql 3> m 2 , 
where Qo is the starting scale for QCD evolution. To find the BFKL asymp¬ 
totics we have to find a different limit of the perturbative series, namely, 

OO 

x b G(x b ,Q' 2 ) = lim| Xfl _»o E a s(Q 2 ) x M n (x B ,Q 2 ) = ( 16 ) 

n =0 

oo 

E a s(Q 2 ) x Um\x B ^o M n (x B ,Q 2 ) • 

n =0 

The change of the order of two operations: sum and limit, looks even more 
suspicious in this case since we have no arguments based on renormalization 
group. For M n we have at x —> 0 

lim\ XB ->o M n (x B ,Q 2 => (17) 

K n (Q 2 ) ln n (l/ x) + K^\Q 2 ) ln n_1 (l/ x) + ... 

Taking the first term in Eq. (17) we arrive to the leading log(l/x) approxima¬ 
tion, namely, 

OO 

x b G(x b ,Q 2 ) = E K n (Q 2 )(a s (Q 2 ) \n(l/x)) n . (18) 

n =0 

Actual summation has been done using the BFKL equation [7] and the answer 
for the series of Eq. (17) is the leading order solution to the BFKL equation 
(LO BFKL). Taking into account term K^(Q 2 ) ln n-1 (l/x) in Eq. (17) we can 
calculate a next order correction to the BFKL equation (NLO BFKL ). 

The LO BFKL asymptotics has a Regge - like form: 



X B G{xBi Q 2 ) 


(19) 



where r = ln(Q 2 /A 2 ) and r 0 = ln(Qo/A 2 ), y = ln(l /xb) and y 0 = 
\w(1/xb,o)- The value of xb,o, which characterize the energy from which we 
can start to use the BFKL asymptotics, cannot be evaluated in the BFKL 
approach. Two parameters ul and D are well defined by the BFKL kernel ( 
see [7] ). 

Two terms in the exponent in Eq. (19) have two different meaning in physics: 
(i) the first one reflects the power - like energy behaviour and, therefore, re¬ 
produces the Reggeon - like exchange at high energy; (ii) while the second 
one describes so called diffusion in log of transverse momentum and, there¬ 
fore, manifest itself a new properties of our microscopic theory - dimensionless 
coupling constant in QCD. 

Below, we will explain in more details all characteristic features of the BFKL 
asymptotics. 

7. Regge factorization: 

Regge factorization ( don’t mix up it with “hard” factorization ) is 

nothing more than Eq. (9) for a Reggeon exchange. Eq. (9) says that the 
amplitude of elastic or quasi-elastic ( like an amplitude for the diffraction 
dissociation reaction ) process can be written as a product 

A(s,t ) = g (projectile, excitation of projectile) (^) (20) 

g ( target, excitation of target ) (££) P ( Reggeon propagator ) . 

In Eq. (20) only Reggeon propagator depends on energy s ( and momentum 
transfer t ), while all dependence on quantum numbers and masses as well as 
on other characteristics of produced particles are factorized out in the product 
of two vertices. It is easy to see that Eq. (20) leads to a large number of 
different relations between measured cross sections. 

I give several examples of such relations to demonstrate how Regge factoriza¬ 
tion works. 

1. the first relation which was suggested by Gribov and Pomeranchuk: 

t) x a tot (pp) = CTt ot (np) . ( 21 ) 


Of course, it is impossible to check Eq. (21) since we cannot measure cr tot (irn) 
but it can be used to estimate its value. 

2 . 


a SD (p+p 


p —> M + p ) B sd 
a el (pp ) B el 


= Const(s) , 


where B is the slope of the cross section in the parameterization 

by _ e -Biii 

f (t = o) 



3. 


a SD (p + p —> M + p) B SD (p + p —> M + p) o el (pp) B el (pp) 

a SD (p + it —> M + 7r) B SD (p + 7r —> M + 7r) cr el (7rp) B el (r:p) 

I think, that these examples are enough to get the spirit what the Regge 
factorization can do for you. You can easily enlarge the number of predictions 
using Eq. (20). 

8. Factorization theorem ( “hard” factorization): 

Any calculation of “hard” processes is based on the factorization theorem 
[8], which allows us to separate the nonperturbative contribution from large 
distances (parton densities, F\(p 2 )) from the perturbative one (“hard” cross 
section, a hard ). For example, the cross section of the high p t jet production in 
hadron - hadron collisions can be written schematically in the form: 

a (A + B —> jets(p t ) + X ) = 

F\(p 2 ) (S) ^A(h 2 ) (7 hard (pa/rtons with p t > k t > /j) . (24) 

cjhard j n ^ 24 ) should be calculated in the framework of perturbative QCD 
and not only in the leading order of pQCD but also in high orders so as 
to specify the accuracy of our calculation. Practically, we need to calculate 
the “hard” cross section at least in the next to leading order, to reduce the 
scale dependence, which appears in the leading order calculation, as a clear 
indication of the low accuracy of our calculation. Of course, we have to adjust 
the accuracy in the calculation of F(p 2 ) and <j hard . 

9. Secondary Reggeons ( trajectories ): 

We call secondary Reggeons or secondary trajectories all Reggeons with inter¬ 
cepts smaller than 1. There is a plenty of different Reggeons with a variety 
of the different quantum numbers. However, they have several features in 
common: 

1. A Reggeon describes the family of resonances that lies on the same Reggeon 
trajectory a^(t). Note, that the Pomeron does not describe any family of 
resonances. In Fig. 2 you can see the experimental p and / trajectories. 

2. All secondary trajectories can be parameterized as 

a R (t) = a R ( 0) + c4(0)t, (25) 

both in scattering ( t < 0 ) and resonance ( t > 0) regions (see Ref. [9] for 
details and for discussion of corrections to Eq. (25) ). It is interesting to note 
that the value of the slope ( cd n (0) ) is the same for all secondary Reggeons 
with the same quark contents . For the Reggeons that corresponds to the 






FIGURE 3. Duality between the resonances in the s - channel and Reggeon exchange in the t - 
channel. 

resonances made of the light quarks a' R (0 = 1 GeV~ 2 ( see Fig. 2). However, 
for Reggeons with heavy quark contents the slope is smaller [9]. 

3. For Reggeons we have duality between the Reggeon exchange and the 
resonance contributions [10] shown in Fig. 3. 

The direct consequence of the duality approach is an idea that any Reggeon 
can be viewed as the exchange of quark - antiquark pair in the t - channel ( see 
Fig.4 ). Therefore, the dynamics of the secondary Reggeons is closely related 
to non-perturbative QCD in quark - antiquark sector while the Pomeron is 
the result of the non-perturbative QCD interaction in the gluon sector. 

10. Shadowing Corrections ( SC ): 

To understand what is the shadowing correction ( SC ) and why we have to 
deal with them let us look at Fig.l more carefully. Actually, the simple formula 
of Eq. (1) is correct only if the flux of “wee” partons N is rather small. Better 
to say, that this flux is such that only one or less then one of “wee” partons has 
the same longitudinal and transverse momenta. However, if N > 1, several 
“wee” partons can interact with the target. For example, two “wee” partons 
interact with the target with a cross section 

cr^ = cr 0 X K (26) 

where k is the probability for the second “wee” parton to meet and interact 
with the target. All “wee” partons are distributed in an area in the transverse 
plane which is equal nR 2 (s) and R 2 (s ) —> a' P (0) ln(s). Therefore, a ( ~ 2> is 
equal to 
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M + B -> M + B 


FIGURE 4. Quark diagrams for meson - meson and meson - baryon scattering. 

= «T 0 X N X K = . (27) 

7 r R z (s) 

Now, let us ask ourselves what should be the sign for a contributions to 
the total cross section. The total cross section is just the probability that an 
incoming particle has at least one interaction with the target. However, in 
Eq. (1) we overestimate the value of the total cross section since we assumed 
that every parton out of the total number of “wee” partons N is able to interact 
with the target. Actually, in the case of two parton interaction, the second 
parton cannot interact with the target if it is just behind the first one. The 
probability to find the second parton just behind the first one is equal k which 
we have estimated. Therefore, instead of flux N of “wee” partons in Eq. (1) 
we should substitute the renormalized flux, namely 

RENORMALIZED FLUX = (28) 

N x (1 - K) = JVx ( 1- °AA- 

V ' V * R2 (s 

The difference between flux of “wee” partons and the renormalized flux of 
“wee” partons if this difference is not very large we call shadowing and/or 
screening corrections. In other words, 

the shadowing corrections for the Pomeron exchange is the Glauber 
screening for the flux of “wee” partons. 

The analogy and terminology become more clear if you notice that Eq. (28) 
leads to Glauber formula for the total cross section, namely, 














(29) 
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a tot 


tot . 


7 r K 2 (s) 


If k <C 1 we can restrict ourselves to the calculation of interaction of two 
“wee” partons with the target. However, if n ps 1 we face a complicated and 
challenging problem of calculating all SC. We are far away from any theoretical 
solution of this problem especially in the “soft” interaction. What we know 
will be briefly reviewed in this paper. 

11. Saturation of the parton ( gluon ) density: 

Saturation of the parton (gluon ) density [11] is a hypothesis about the renor¬ 
malized flux of “wee” partons in the region where flux N is much large than 
unity. More precisely, we assume that the value of n reaches the unitarity 
maximum ( k — 1) and ceases to increase. Fig. 5 gives a picture of the par- 
ton density in a target which corresponds this hypothesis. From this picture 
one sees that the unitarity limit can be reached at sufficiently large virtualities 
Q 2 ( short distances ) in the region of applicability of pQCD. It allows us to 
evaluate better the flux N and parameter k. 

The expression for k can be written in the form 

k = xG M ^ = xG(x, Q 2 ) , (30) 


where 

1. xG(x,Q 2 ) = the number of partons ( gluons) in the parton cascade (IV); 

2. R 2 is the radius of the area populated by gluons in a nucleon; 

3. a(GG) is the gluon cross section inside the parton cascade and was evalu¬ 
ated in Ref. [12], 

Actually, at present we know xG and R 2 well enough [44] to estimate the 
kinematic region where the flux of “wee” partons becomes large. This is the 
line (see Fig.5 ) 

n(x,Q 2 ) = 1. (31) 

We will discuss below the physics of this equation. 

12. Semiclassical gluon field approach: 

In the region where iV > 1 we can try to develop a different approach 
[13]- it semiclassical gluon field approach. Indeed, due to the uncertainty 
relation between the number of particles and the phase of the amplitude of 
their production 


N x Af > 1 , 


(32) 
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FIGURE 5. Parton distributions in the transverse plane as function of ln(l/x) and Q 2 . 
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FIGURE 6. Impact parameter b t and useful kinematic relations for high energy scattering. 


we can consider Aef) oc | < 1. Therefore, we can approach such parton 
system semiclassically. It means that the parton language is not more suited 
to discuss physics and we have to find an effective Lagrangian formulated in 
term of semiclassical gluonic field. A remarkable theoretical progress has been 
achieved in framework of such an approach [13] but we mention this rather 
theoretical approach here only to illlustrate the richness of ideas that brought 
to the market the experiments in the deep inelastic processes mostly done at 
HERA and the Tevatron. 

13. Impact parameter ( b t ) representation: 

For high energy scattering it turns out to be uscfull to introduce impact pa¬ 
rameter b t by 

l = k x b t , 

where l is the angular momentum in s - channel of our process and k is the 
momentum of incoming particle at high energy. High energy means that we 
consider the scattering process for fca > 1 where a is the typical size of the 
interaction. 

The meaning of b t as well as useful kinematic relations for high energy scat¬ 
tering is pictured in Fig.6. 

Formally, the scattering amplitude in 6 t -space is defined as 



a e i(s,b t ) = J d 2 q e tql ' bt A(s,t) 

(33) 

where A(s,t) is scattering amplitude and t = —q\ is 
squared. The inverse transform for A(s,t) is 

momentum transfer 

A(s,t) = J aei(s, b t ) d\ e iq± ' bt . 

(34) 

In this representation 


o lol = 2jd%Ima el (sM; 

(35) 

a d = J d\ \a e i(s,b t )\ 2 . 

(36) 


14. s - channel unitarity: 

The most important property of the impact parameter representation is the 
fact that the unitarity constraint looks simple in this representation, namely, 
they can be written at fixed b t . 

2 Inn a e i(s, b) = \a d (s,b)\ 2 + G in (s,b) , (37) 

where G m is the sum of all inelastic channels. This equation is our master 
equation which we will use frequently during this talk. 

15. A general solution to the unitarity constraint: 

Our master equation ( see Eq. (37) ) has the general solution 

G in (s,b) = 1 - e- n( ^ ; (38) 

a el = i{l - g-^+ixWO} ; 

where opacity 0 and phase y = 2 5(s, b ) are real arbitrary but real functions. 

The algebraic operations, which help to check that Eq. (36) gives the solution, 
are: 

_ n 

Ima e i = 1 — e 2 cosy ; 

|a e ;| 2 = (l — e“§ + lx )-(l — e~^“* x ) = 1 — 2e“§ cosy + e~ n . 

The opacity Q has a clear meaning in physics, namely e~ n is the probability 
to have no inelastic interactions with the target. 

One can check that Eq. (37) in the limit, when O is small and the inelastic 
processes can be neglected, describes the well known solution for the clastic 



scattering: phase analysis. For high energies the most reasonable assumption 
just opposite, namely, the real part of the elastic amplitude to be very small. 
It means that X —^ 0 and the general solution is of the form: 

G m (s,b) = 1 - e- n( ^ ; (39) 

_ n(a,b) 

CL e l = l { 1 - 6 ^ } . 

We will use this solution to the end of this talk. At the moment, I do not 
want to discuss the theoretical reason why the real part should be small at 
high energy . I prefer to claim that this is a general feature of all experimental 
data at high energy. 

16. The great theorems: 

Optical theorem 

The optical theorem gives us the relationship between the behaviour of the 
imaginary part of the scattering amplitude at zero scattering angle and the 
total cross section that can be measured experimentally. It follows directly 
from Eq. (37), after integration over b. Indeed, 

4nImA(s, t = 0) = J 2Ima e i(s,b t ) d 2 b t = (40) 


J d 2 bt{ \a e i(s, 6j) | 2 + Gi n (s,bt ) — cr cl + cr in — cr tot ■ 

Froissart bound 

We call the Froissart boundary the following limit of the energy growth of the 
total cross section: 


(Jtot < C In 2 s (41) 

where s is the total c.m. energy squared of our clastic reaction: a{p a ) +b(pb) —> 
a + b , namely s = (p a + pb) 2 . The coefficient C has been calculated but we 
do not need to know its exact value. What is really important is the fact 
that C oc A where kt is the minimum transverse momentum for the reaction 
under study. Since the minimum mass in the hadron spectrum is the pion 
mass the Froissart theorem predicts that C oc -K. The exact calculation 
gives C = 60 mb. 

The proof [14] is based on Eq. (39) and on the proven asymptotics for 12, 
namely, 


bt) Ib t > i 


s N e - 2 » bt 


(42) 


where /j is the mass of the lightest hadron ( pion). 



It consists of two steps: 

1. We estimate the value of 6° from the condition 

SIM?) « 1 • (43) 


Using Eq. (42) one contains 


6° = — Ins 

* 2 fi 


(44) 


Note, that at high energies the value of 6? does not depend on the exact value 
of n in Eq. (43). 

2. Eq. (35) for the total cross section we integrate over b t by dividing the 
integral into two parts b t > and b t < 6°. Neglecting the second part of the 
integral where is very small yields 



This is the Froissart bound. 

Pomeranchuk theorem: 

The Pomeranchuk theorem is the manifestation of the crossing symmetry, 
which can be formulated in the following form: one analytic function of two 
variables s and t describes the scattering amplitude of two different reactions 
a + b —> a + & at s > 0 and t < 0 as well as a + b —> a + b at s < 0 
(u = (p a + Pb ) 2 > 0 ) and t < 0. 

The Pomeranchuk theorem says that the total cross sections of the above two 
reactions should be equal to each other at high energy 

&tot(a + b) = a tot (a + b) at s -> oo , (45) 


if the real part of the amplitude is smaller than imaginary part. 

17. The “ black” disc approach: 

The rough realization of what we expect in the Froissart limit at ultra high 
energies is so called “black disc” model, in which we assume that 


= oo at b < R(s ) ; 
= 0 at b > R(s ) . 



The radius R can be a function of energy and it can rise as R ~ In s dne to 
the Froissart theorem. 

It is easy to see that the general solution of the unitarity constraints simplifies 
to 

a(s, b) e i = i®(R — b) * 

G in (s,b ) = Q(R - b) ; 

which leads to 


(Tin f d bC in 

= 7 T R 2 ] 

(46) 

a et = J d 2 b\a e i\ 2 

= n R 2 ; 

(47) 

CTtot = a el + Gin 

= 2ir R 2 ; 

(48) 


A(s,t) = £ / J 2 bO(R - 6 )e‘* s = itf bdbMbq,) = i§Mq,R) ; 

§ = *\A \ 2 = , ( 49) 

Comparison with the experimental data of the “black disc” model shows that 
this rather crnel model predicts all qualitative features of the experimental 
data such as the value of the total cross section and the minimum at certain 
value of t, furthermore the errors in the numerical evaluation is only 200 - 
300%. This model certainly is not good but it is not so bad as it could be. 

Actually, this model as well as the Pomeron approach is the first model that 
we use to understand the new experimental data at high energy. 

18. Generalized VDM for photon interaction: 

Gribov was the first one to observe [15] that a photon ( even virtual one) 
fluctuates into a hadron system with life time ( coherence length ) r = l c — 
where Xb = Q 2 is the photon virtuality and m is the mass of the 
target. This life time is much larger at high energy than the size of the target 
and therefore, we can consider the photon - proton interaction as a processes 
which proceeds in two stages ( see Fig.7): 

(a) Transition 7 * —> hadrons which is not affected by the target and, there¬ 
fore, looks similar to electron - positron annihilation; 

(b) hadron — target interaction, which can be treated as standard hadron 
- hadron interaction, for example, in the Pomeron ( Reggeon ) exchange 
approach . 

These two separate stages of the photon - hadron interaction allow us to use 
a dispersion relation with respect to the masses M and M' [15] to describe 
the photon - hadron interaction ( see Fig .8 for notations ), as the correlation 
length l c = A > Rv, the target size. Based on this idea we can write a 
general formula for the photon - hadron interaction, 
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FIGURE 7 . Two stages of photon - hadron interaction at high energy. 
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(50) 


where T 2 (M 2 ) = R(M 2 ) (see Eq. (52) ) and cr(M 2 ,M /2 ,s) is proportional to 
the imaginary part of the forward amplitude for V + p —» V' + p where 

V and V are the vector states with masses M and M'. For the case of 
the diagonal transition ( M = Ad') cr(M 2 ,s) is the total cross section for 

V — p interaction. Experimentally, it is known that a diagonal coupling of 
the Pomeron is stronger than an off-diagonal coupling. Therefore, in Erst 
approximation we can neglect the off-diagonal transition and substitute in 
Eq. (50) cr(M 2 , M /2 , s) = a)M 2 .s)M 2 5(M 2 - Ad' 2 ). 

The resulting photon - nucleon cross section can be written written as: 


cr(7*iV) 


a em f R(Ad 2 )M 2 dM 2 

J (Q 2 + Ad 2 ) 2 


(51) 


where R(Ad 2 ) is defined as the ratio 


R(Ad 2 ) 


cr(e + e —> hadrons) 
cr(e + e~ —> p + p~) 


(52) 


The notation is illustrated in Fig.8 where Ad 2 is the mass squared of the 
hadronic system, T 2 (M 2 ) = R(Ad 2 ) and a M 2 N (s) is the cross section for the 
hadronic system to scatter off the nucleonic target. 

Experimentally, R(M 2 ) can be described in a two component picture: the con¬ 
tribution of resonances such as p,u, 0,J/T, T' and so on and the contributions 
from quarks which give a more or less constant term changes abruptly with 
every new open quark - antiquark channel, R(Ad 2 ) « 3 J2 q e q i where e q is the 
charge of the quark and the summation is done over all active quark team. 




FIGURE 8. The generalized Gribov’s formula for DIS. 


If we take into account only the contribution of the p - meson, lu - meson , <f> 
and J/'k resonances in R(M 2 ) we obtain the so called vector dominance model 
( VDM) [16] which gives for the total cross section the following formula: 


Clot (7* + P) 


a P 


3n 


vtotip + P)J2 R ( M2 = m u){ 




Q 2 + < 


(53) 


where my ; is the mass of vector meson, Q 2 is the value of the photon virtuality 
and i?(M 2 = rriy) is the value of the R in the 

mass of the vector meson which can be rewritten through the ratio of the 
electron - positron decay width to the total width of the resonance. Of course, 
the summation in Eq. (53) can be extended to all vector resonances [17] or 
even we can return to a general approach of Eq. (50) and write a model for 
the off-diagonal transition between the vector meson resonances with different 
masses [18]. We have two problems with all generalization of the VDM: (i) a 
number of unknown experimentally values such as masses and electromagnetic 
width of the vector resonances with higher masses than those that have been 
included in VDM and (ii) a lack of theoretical constraints on all mentioned 
observables. These two reasons give so much freedom in fitting of the experi¬ 
mental data on photon - hadron interaction that it becomes uninteresting and 
non - instructive. One can see from Eq. (53), that if Q 2 my VDM predicts 
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FIGURE 9-a. Lego - plot for double diffractive dissociation. 

a ^ behaviour of the total cross section. Such a behaviour is in clear contra¬ 
diction with the experimental data which show an approximate ^ dependence 
at large values of Q' 2 ( i.e. scaling ). 

The solution of this puzzle as well as the systematic description of the pho¬ 
ton - proton interaction can be reached on the basis of Gribov formula but 
developing a general description both “soft” and “hard” processes. The first 
and oversimplified version of such description was suggested by Bodelek and 
Kwiecinski [19] and more elaborated approaches have been recently developed 
by Gotsman, Levin and Maor [20] and Martin, R.yskin and Stasto [21], 

19. Difftraction dissociation: 

The experimental definition of the diffractive dissociation processes: 

Diffractive dissociation processes are processes of production one ( 
single diffraction (SD)) or two groups of hadrons ( double diffraction 
( DD ) ) with masses ( and M 2 in Fig.9 ) much less than the 
total energy ( Mi <C s and M 2 <C s ). 

In diffractive processes no hadrons are produced in the central region of ra¬ 
pidity as it shown in Fig.8. This is the reason why diffractive dissociation is 
the simplest process with large rapidity gap ( LRG ). 

The above definition of diffractive processes is the practical and/or experimen¬ 
tal one. From theoretical point of view as was suggested by Feinberg [22] and 
Good and Walker [23] diffractive dissociation is a typical quantum mechanical 
process which stems from the fact that the hadron states are not diagonal with 
respect to the strong interaction. Let us consider this point in more details, 
denoting the wave functions which are diagonal with respect to the strong in- 



teraction by T n . Therefore, the amplitude of high energy interaction is equal 
to 

A n = < V n \T\y n , > = A n 8 ny , (54) 

where brackets denote all needed integration and T is the scattering matrix. 
The wave function of a hadron is equal to 

OO 

hadron = • (55) 

?i=l 

Therefore, after collision the scattering matrix T will give you a new wave 
function, namely 

OO OO OO 

^ final E £Cn<*„ |T|T„/ > T n , = Y. C nAn^n- (56) 

n =1 n '=1 n= 1 

One can see that Eq. (56) leads to elastic amplitude 

OO 

O'el = ^ ^ final \ ^hadron = ^ ] C n A n , ( 57 ) 

n =1 

and to another process, namely, to production of other hadron state since 
T final ^ ^hadron- This process we call diffractive dissociation. The total 
cross section of such diffractive process we can find from Eq. (56) and it is 
equal to 

O (s, bt) = 'y ] < final \ ^hadron ^ ^hadron \ ^hadron = (58) 

final 

oo / oo \ 2 

E Cl Alls, b t ) - E C?A.(»,6i) , 

n=l \ n=l / 

where we regenerate our usual variable: energy (s) and impact parameter ( b t ). 
Using the normalization condition for the hadron wave function ( E n = 1 
) we can see that Eq. (58) can be reduced to the form [23] [24] 

o SD (s,b t ) = < \a 2 (s,b t )\ > - < \o(s,b t )\ > 2 , (59) 

where <[/[>= E n <E./n- 

Eq. (59) has clear optical analogy which clarify the physics of diffraction, 
namely, the single- slit diffraction of “white” light. Indeed, everybody knows 
that in the central point all rays with different wave lengths arrive with the 
same phase and they give a “white ” maximum. This maximum is our elastic 
scattering. However, the conditions of all other maxima on the screen with 
displacement a, depend on the wave lengths ( asind — k A ). All of them are 
different from the central one and they have different colours. Sum of all these 
maxima gives the cross section of the diffractive dissociation. 



20. Pumplin bound for SD: 

The Pumplin bound [24] follows directly from unitarity constrain of Eq. (37) 
which we can write for each state n with wave function T r) separately: 

2 ImAi(s,b t ) = K‘M,)| 2 + ■ (60) 

Assuming that the amplitude at high energy is predominantly imagine, we 
obtain, as has been discussed, that 

| 2 < G™(s,b t < a^(s, b t ) . (61) 

After summing over all n with weight C 2 ( averaging ) one obtains the Pumplin 
bound: 


Vei(s,b t ) + s SD (s,b t ) < a tot (s,b t ) , (62) 

where a e i(s.b t ) = \a e i(s,b t )\ 2 , a to t(s,b t ) = 2 Im a e ;(s, b t ) and s SD (s, b t ) = 
Em | a M (s,b t ) f . 

21. Unitarity limit of SD for hadron and photon induced reactions: 

At high energy A ==>• 1 ( see Eq. (41) ). It means that at high energy 
the scattering matrix does not depend on n. Therefore from Eq. (56) one can 
derive 


-final 


(63) 


y T. °n < > i i>„/ = y c„4i„ = . 

n= 1 n '=1 n= 1 

Eq. (63) says that cr SD 0 at high energies. Performing integration over 

bt it is easy to obtain that in the hadron induced reactions 


<7 


SD 


&el 



(64) 


It should be stressed, that we should be careful with this statement for photon 
( real or virtual ) induced reactions. As we have seen, in Gribov formula we 
have an integration over mass. For each mass we have the same property of 
the diffractive dissociation as has been presented in Eq. (64). However, due to 
integration over mass in Eq. (51)the ratio of the total diffractive dissociation 
process to the total cross section is equal to 


-> M + P) 
<?tot(l*P -»■ 7 * +P) 


1 

2 


(65) 



22. Survival Probability of Large Rapidity Gaps: 

We call any process a large rapidity gap ( LRG )process if in a sufficiently large 
rapidity region no hadrons are produced. Historically, both Dokshitzer et al. 
[25] and Bjorken [26], suggested LRG as a signature for Higgs production in W- 
W fusion process in hadron-hadron collisions. The definition of the survival 
probability ( < S 2 > ) is clear from Fig. 9-b. Indeed, let us consider the 
typical LRG process - production of two jets with large transverse momenta 
Pti ~ — Pt .2 ^ G where fi is the typical mass scale of “soft” processes , and 
with LRG between these two jets. Therefore, we consider the reaction: 

p( 1 ) + p( 2 ) —» M^hadrons + jefi(yi,Pti)] ( 66 ) 

+ LRG[Ay = |j/i - y 2 |] + M 2 [hadrons + jef 2 (y 2 ,pt 2 )] , 

where y\ and y 2 are rapidity of jets and LRG Ay — \y 1 — y 2 \ 1. To produce 

two jets with LRG between them we have two possibility: 

(a) This LRG appears as a fluctuation in the typical inelastic event. However, 

Ay 

the probability of such a fluctuation is proportional to e~~ and the value 
of the correlation length L we can evaluate L ~ where ^ is the 

number of particles per unit in rapidity. Therefore, LRG means that 

Ay » L ; 

(b) The exchange of colourless gluonic state in QCD is responsible for the 
LRG. This exchange we denote as a Pomeron in Fig.9-b. The ratio the 
cross section due to the Pomeron exchange to the typical inelastic event 
generated by the gluon exchange ( see Fig. 9-b ) we denote F s . Using the 
simple QCD model for the Pomeron exchange, namely, the Low-Nussinov 
[27] idea that the Pomeron = two gluon exchange, Bjorken gave the first 
estimate for F s « 0.15%. This is an interesting problem to obtain a 
more consistent theoretical estimates for F s , but F s is not the survival 
probability. F s is the probability of the LRG process in a single parton 
shower. 

However, each parton in the parton cascade of a fast hadron could create its 
own parton chain which could interact with one of the parton of the target. 
Therefore, we expect a large contribution of inelastic processes which can fill 
the LRG. To take this multi shower interaction we introduce the survival 
probability < S 2 > ( see Fig.9-b ). 

To calculate < S 2 > we need to find the probability that all partons with 
rapidity > yi in the first hadron ( see Eq. (66) and partons with yj < y 2 
in the second hadron do not interact inclastically and, therefore, they cannot 
produce hadrons in the LRG. As we have discuss, the meaning of Eq. (41) in 
physics is that 

P(yi- yj M j ) = e- n <»*-w^) (67) 



gives you the probability that two partons with rapidities y % and y 3 and with 
b\ 3 = | b t) i — b t j | do not interact inelastically. 

The general formula for < S 2 > reads 

S 2 / dyj dyj d 2 b\ 3 Pjjy^ y u 6 M ) P 2 (yj, y 2 , b t J) Pijh ~ Vj, K 3 ) , . 

f dyi dy 3 d 2 b\ 3 P 1 (y t , y x ) P 2 (yj, 2 / 2 ) 

where P\ (y ? , y\){P-ziy^ IJ 2 ) is probability to find two partons with rapidity y\ 
and yi ( 7/2 and y 3 ) in hadron 1(2), respectively. The deep inelastic structure 
function is x 1 G(x 1} p 2 1 ) = J dyiCpb t>i Pi(y 1 = ln(l/xi), y u b t>i ). 

Unfortunately, there is no calculation using a general formula of Eq. (68). 
What we have on the market is the oversimplified calculation in the Eikonal 
model ( see Ref. [28] ) in which we assume that only the fastest partons from 
both hadrons can interact. 
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FIGURE 9-b. Typical LRG process . 
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FIGURE 10. The exchange of the Pomeron. 

I “S O F T” POMERON: 

II “THEOREMS” ABOUT POMERON 

In this section we collect main properties of the “soft” Pomeron which we have 
discussed in the previous section. 

• There is only one Pomeron. 

It means that only “soft” Pomeron has a chance to be a simple Regge pole 
(Reggeon) ( see previous section for details). 

• The Pomeron is the Reggeon with A P (0) = a P (0) — 1 <C 1 

Being a Reggeon, the Pomeron has all typical features of a Reggeon written 
in Eq. (9) and presented in Fig. 10. They are: 

-k Analiticity ; 

-k Resonances ; 

* Factorization ; 

* Definite phase ; 

-k Increase of the radius of interaction ; 

k Shrinkage of the diffraction peak ; 

Increase of the radius of interaction: 

Eq. (9) can be rewritten in the impact parameter representation taking the 
integral of Eq. (33) as follows 



(69) 


a e i(s,b t ) = 9i(0)g;(Q)(s/so) Ap ^7^y e ”’<•> 

with 

i? 2 (s) = i?g + 4 ap ln(s/s 0 ) 

One can see, that the typical impact parameters turns out to be of the order 
of R(s). This fact means that the radius of interaction increases with energy. 
Shrinkage of the diffraction peak: 

Using Eq. (9) we can see that the f-dependence is determined mostly by the 
Pomeron exchange since 


da 

dt ^ e -2a' p \n{s/s 0 )\t\ . 

W l*=« 

Therefore, the differential cross section falls down rapidly for \t\ > 

a' (o) in(s/s 0 ) • This phenomena we call the shrinkage of the diffraction peak. 
It should be stressed that this shrinkage has been observed in many “soft” 
reactions, but it seems to me that the value of a' P (0) has not been determined 
yet from the experimental data. 

Pomeron = gluedynamics at high energy . 

There is no resonances on the Pomeron trajectory. This is the principle differ¬ 
ence between the secondary Reggeons and the Pomeron. In duality approach 
the Pomeron does not appear in dual diagrams of the first order and 


P oc 


as_ 

N c 


S \Ca s 
s 0 


The common believe, which is 

based on the experience with duality approach and QCD calculations, is that 
The Pomeron —> glueballs —» gluon interaction at high energy. 


There is no arguments for a Pomeron but it describes experimental 
data quite well 

Fig. 11 shows that the Pomeron exchange is able to describe the experimental 
data on total proton - proton scattering at high energy. Of course, this is 
an example. More systematic description of the data one can find in Ref. [4] 
or/and in original Donnachie - Landshoff papers [3]. 

Pomeron is a “ladder” diagram for a super convergent theory like 
g 0 3 - 

The charm of the Pomeron approach is in the relation between the elastic 
or / and quasi - elastic reactions ( like diffractive dissociation, for example ) 
and the multiparticle production at high energy. The simple picture for the 




FIGURE 11. The total cross section in the Pomeron approach in the Donnachie - Landshoff 
approach (full line ) and in the Eikonal model (dotted line) 
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FIGURE 12. The Pomeron structure in gcf> 3 - theory ( parton model ). 


Pomeron exchange presented in Fig. 12 was, is and, unfortunately, will be our 
practical tool for application of the Pomeron phenomenology to the processes 
of the multiparticle production. 

In spite of the simplicity of this picture we would like to mention that it appears 

in simple but strict approximation, so called leading log s approximation. In 

this approximation we sum all contributions to Feyman diagrams of the order 
2 

of ( ^2 hi s ) n , neglecting smaller terms. 

It should be stressed that this approach reproduces the main property of ob¬ 
served multiparticle production at high energies. 

-k Power - like behaviour of the total cross section 

a tot oc (-f^ ; 
so 


-k Pomeron = only multiparticle production with < N > oc In .s ; 

* Uniform rapidity distribution of the produced particles ; 

-k Average transverse momentum of produced particles does not depend on 
energy ; 

-k There are no correlation between produced particles ; 
k Poisson multiplicity distribution ; 

* Increase of the interaction radius due to the random walk in impact param¬ 
eters ; 
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FIGURE 13. Uniform distribution in rapidity in g<jr - theory ( parton model ) . 



FIGURE 14. Random walk in the transverse plane in g(f> 3 - theory ( parton model ). 
























The increase of the interaction radius is seen in Fig. 14 and follows directly from 
uncertainty relation and from the fact that the mean transverse momentum of 
particles ( partons ) does not depend on energy. Indeed, from uncertainty principle: 

A b t x < kt > « 1 , (71) 


where < kt > is the mean parton momentum. Therefore, each emission changes 
the position of the parton in impact parameter space on the value A bt ~ <j ^ . 
After n emission the parton will be on the distance < > oc < ^ > n. Since the 

average number of emission N oc In s, the radius of interaction 


R 2 (s ) = < bf N > oc 


< k t > 


■N = 


<k f > 


In s . 


(72) 


III DONNACHIE - LANDSHOFF POMERON 

In practice, when we say “soft” Pomeron, we mean so called Donnachie - Land- 
shoff Pomeron. Donnachie and Landshoff suggested the simplest picture of high 
energy interaction - the exchange of the Pomeron which is a Regge pole. Note, 
that the Donnachie - Landshoff approach to high energy scattering is 
more complicated than the exchange of the Pomeron. They also in¬ 
cluded a week shadowing corrections as we will discuss below. I want to 
stipulate that in this section I discuss only the D-L Pomeron but not the D-L ap¬ 
proach. In their simple picture Donnachie and Landshoff gave the most economic 
and elegant description of the experimental data mostly related to the total and 
clastic cross sections. From fitting procedure they found that the D-L Pomeron has 
the following features: 

1. The Pomeron intercept, A = ctp(O) — 1 = 0.08 ; 

2. The Pomeron slope ap(0) = 0.25 GeV~ 2 ] 

3. The linear parameterization for the Pomeron trajectory otp(t) = CKp(O) + 
CKp(0) t can be used to fit the experimental data ; 

4. The additive quark model ( AQM ) can be used to to find the t - dependence 
of the Pomeron - hadron vertices. In AQM the Pomeron vertex is proportional to 
the electromagnetic form factor of the hadron. 

IV SEVEN ARGUMENTS AGAINST D - L POMERON 

In spite of the wide use of the D-L Pomeron or, may be, because of this, the 
D-L Pomeron is seriously sick. Here, we want to list seven arguments against the 
D-L Pomeron which, we think, show that D-L Pomeron ( and D-L approach, which 
we will discuss below) is close to its last days of existence. In what follows we 
try to separate the D-L Pomeron from D-L approach which includes not only the 
Pomeron but also week SC. 
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FIGURE 15. a e i for the D-L Pomeron . 


1. The D-L Pomeron violates unitarity ( a e i (s,b t ) < 1 just around the 
corner ) ; 

From Fig.15 one can see that at yfs ~ 3000GeP the D-L Pomeron violates the 
s-channel unitarity constrain ( see Eq. (37) ). The success of the description 
of the experimental data on total and elastic cross section is mostly due to the 
fact that the area, in which a e i(s, bt) > 1 in Fig. 15, is much smaller than the 
total area. However, the fact that a e i is close to unity should reveal itself in 
the description of other processes which are more sensitive to the value of a e i. 

2. The D-L Pomeron gives Gi n (s, b t ) which is quite different from the 
Pomeron one ; 

In Fig. 16 we plot the value of Gi n (s, b t ) which we calculated from the s- 
channel unitarity constrain of Eq. (37), namely, 


G in (s,b t ) = 2 Im a et (s, b t ) - \a e i(s.b t )\ 2 , 

using the D-L parameterization for a e i(s,b t ). One can see that Gi n turns out 
to be very close to unity in accessible energy range. Comparing ^-dependence 
of a e i(s,b t ) ( see Fig.14 ) and G in (s,b t ) ( see Fig.15) one can see that this 
dependence is quite different for them. It means that the main idea of the 
parton model (see Fig.l) that the exchange of the Pomeron is closely related 
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FIGURE 16. Gi n for the D-L Pomeron . 


to the multiparticle ( multiparton ) production is broken in the D-L approach. 
In other words, the most charming feature of the Pomeron approach as a 
whole, namely, the ability of the Pomeron to describe both elastic and inelastic 
processes, turns out to be inconsistent in the D-L approach. For me, it is too 
high price for the fit of the clastic and total cross sections even if it is a simple 
one. 


3. The D-L Pomeron and D-L approach cannot describe the s- 
dependence of the diffractive dissociation cross section, measured 
at CDF ; 

For D-L Pomeron a SD oc () 2Ap while a tot oc ( ) Ap . Therefore, for the 
D-L Pomeron 


s 


SD 


&tot 


OC 



(73) 


while experimentally [29] the energy dependence of this ratio is quite different 
as Fig. 17 demonstrates. Week SC of the D-L approach cannot change this 
conclusion. It turns out that SC should be as strong as in the Eikonal model 
or even stronger to describe the CDF data on SD. 

4. The parameters of the D-L Pomeron do not fit data quite well ; 







The parameters of the D-L Pomeron were mostly fitted from the scattering 
data but the D-L Pomeron trajectory predicts a resonance at m 2 ~ 4 GeV 2 
with spin 2 which was not found experimentally. On the other hand, the 
experimental data, especially new HERA data on photon - proton scattering 
show that the shrinkage of the diffraction peak ( «/(()) ) is different for different 
particles [30]. 

5. The D-L Pomeron and /or D-L Approach cannot describe the CDF 
data on double parton cross section ; 

We will comment on this below after discussion of the shadowing corrections. 
We will show that the double parton cross section is closely related to the 
collision of the two parton showers which are absent in the D-L Pomeron and 
too small in the D-L approach. 

6. The D-L Pomeron and /or D-L Approach predicts Survival Prob¬ 
ability for LRG processes —> 1 and, therefore, cannot describe the 
measured LRG processes at the Tevatron ; 

As we have discussed in section 1 the origin of the survival probability is the 
fact that there are many parton shower interactions with the target. If we have 
only the exchange of the Pomeron or, in other words, only;y one parton shower 
interaction the survival probability is equal to unity. Experimentally [31], this 
survival probability is, at least, 0.1 and it shows a substantial s - dependence 
( see Ref [28] for details ). In D-L approach this survival probability is close 
to unity since the value of the SC that they suggested is too small. 


7. The D-L Pomeron and /or D-L Approach cannot reproduce the 
Glauber Shadowing for hadron - nucleus interactions ; 

It is well known that the Glauber Shadowing can be described in the Reggeon 
approach as the multiPomeron exchange of the incoming particle with nu¬ 
cleons in a nucleus. In D-L approach such multiPomeron exchange is sup¬ 
pressed. On the contrary for nucleus - nucleus interaction D-L approach gives 
the Glauber formula since the exchange of the Pomerons between different nu¬ 
cleons from different nuclei are not suppressed. Therefore, in the D-L approach 
the Glauber shadowing for hadron - nucleus and nucleus - nucleus interactions 
look quite different. In the extreme case of the D-L Pomeron approach they 
predict that 

<Jtot(hadron + nucleus) oc A 


while 


a to t(nucleus + nucleus) oc As . 


V SHADOWING CORRECTIONS ( SC ) 


In this section I discuss the shadowing ( screening, absorption ) corrections to 
the Pomeron exchange. Formally speaking,in the Reggeon approach such correc- 



tions can be pictured as the exchange of many Pomerons (see Figs.18a - 18b for 
examples). The only question, that we should answer, is how can we sum all these 
complicated diagrams? 


A Several general remarks: 

• There is no Pomeron without SC ; 

It means that there is no any theoretical idea why many Pomeron exchange 
can be equal to zero. Therefore, our main strategy in operating with the SC is 
to find the kinematic region in which the SC are small, to develop a technique 
how to calculate the SC when they are small and to built a model ( better 
theory of course but practically this is a problem for future) to approach SC 
in the kinematic region when they are large. 

• SC follow from the s-channel unitarity ; 

To understand why SC follow from the s-channel unitarity let us generalize 
the solution to the unitarity constrains of Eq. (37) using the hadronic states 
which are diagonal with respect to strong interaction and which have wave 
functions T n , as we did discussing the diffraction dissociation processes in the 
previous section ( see Eq. (54) - Eq. (59) ). For each of this state we have a 
unitarity constrain of Eq. (??) and we can find a solution to Eq. (60) assuming 
that ImA n (s, b t ) 3> ReA n (s,b t ): 

A eln = i | { 1 — e "V t] } ; (74) 

G in = i _ e -n n (s,b t ) . ( 75 ) 

Let us assume that fl n -C 1 and expand Eq. (74) and Eq. (75) with respect 
to . 

«s,6 ( ) = + O(flJ); (76) 

G“(«A) = Si„(s,b.) - SiifM + O(SJJ). (77) 

Using Eq. (55) - Eq. (58) we obtain for the observables: 



Oei(sM = i (E“iC^„(s,6,)) 2 ; 

(78) 

&t.ot 

= CJO.M,) - iK, Cldl(s,b t )) ; 

(79) 

@in 

= E”i c„ 2 si„M«) - i (E“, claHsM)) ; 

(80) 

SD 

= HE^CnSiJM,) - (ESL, c-=n„( s ,6 t )) 2 }. 

(81) 


As we have discussed, in the framework of converged theories or in the parton 
model, the one Pomeron exchange corresponds to the typical inelastic event 






FIGURE 18-b. The space - time picture in the parton model for the first “fan” Pomeron dia¬ 
gram . 

with production of the large number of particles. Therefore, we can associate 
this exchange with Q since one can see that cr tot (s, = a i n (s,b t ) oc fk All 
terms which are proportional to fl 2 describe the two Pomeron exchange and 
they induce the SC. 

• The scale of SC from the experimental data on Ot ot ,(X e i and cr SD ; 

We can evaluate the scale of the SC using experimental data on a to t, &ei and 
a SD . Indeed, we can write the expression of the total cross section in the form: 

°tot = °tot ~ A<t m > ( 82 ) 

where a[ ot is the contribution of the Pomeron exchange to the total cross 
section. Summing Eq. (78) and Eq. (81) we derive that 

A ^ + ° SD (83) 

or 

_ A a tot _ a e i + a SD f d\ 

crtot o- tot ^ Jd 2 b t n(s,b t ) 

Fig.17 shows that in the wide range of energy n rs 0.34 and, therefore, we 
can consider the single Pomeron exchange as a good approximation only with 
the errors of about 34%. I do not think, we need to argue that such an 
approximation cannot be considered as a good one. 



• The scale of SC from the inclusive correlation function ; 

It is well known ( see for example Ref. [2] ) that the Reggeon approach gives 
the two particle rapidity correlation function, which defined as: 

d 2 cr{yi,y2) 

T3 _ vtotdyidyz 

dajyi) da(y 2 ) 
o'tot dyi o to tdy2 



where 


d 2 a 

dyidy 2 


is the double inclusive cross section for the reaction: 


a + b —> l(?/i) + 2 ( 7 / 2 ) + anything . 


The correlation function can be written in the form: 

At/ 

R(jji, 7/2) = SR x + LR . ( 86 ) 

These two terms in Eq. ( 86 ) correspond to two Reggeon diagrams in Fig.19. 
The first diagram gives the short range correlations which falls at Ay = \y\ — 
7/2 1 > L cor with L cor = a \ Q . i where 0 : 77 ( 0 ) is the intercept of the secondary 
Reggeon trajectory ( see Fig.2 ) and 077 ( 0 ) ~ 0.5. The second term is closely 
related to the shadowing correction to the total cross section and it is equal: 

LR = 2 . (87) 

O'tot 

Therefore, if we can separate the long range rapidity correlation from the short 
range one we will measure the value of the SC. 

The CDF collaboration at the Tevatron did this. The CDF measured the 
process of inclusive production of two “hard” jets with large but almost com¬ 
pensating transverse momenta in each pair( p n cs — p t2 p, where p is the 
scale of “soft” interactions) and with values of rapidity that are very similar. 
Such pairs cannot be produced by one Pomeron exchange or in other words in 
one parton shower collision if the difference in rapidity of these pairs is small 
that l/otsipt)- They can only be produced in double parton shower interac¬ 
tion and their cross section can be calculated using the Mueller diagrams [32] 
given in Fig.20. It can be written in the form [33] 

^ind( 2 jets ) <Tj nc ;( 2 jets ) ,. 

cr D p = m --- , fooj 

2 a eff 

where factor m is equal 2 for different pairs of jets and to 1 for identical pairs. 
The experimental value of cr e ff = 14.5 ± 1.7 ± 2.3 mb [33]. 

Comparing Eq. ( 88 ) with Eq. (87) one can obtain the estimates for 

AtJtot cr m 

Otot 4 (7 e ff 




40 -50% . 


(89) 




FIGURE 20. The Mueller ( Reggeon ) diagram for double parton interaction . 
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FIGURE 21. The first order SC oc J d 2 btQ 2 (s,bt) in the additive quark model for the proton . 

• The scale of SC from diffractive dissociation at HERA ; 

As we have discussed the value of the SC crusially depends on the size of the 

p 

target ( see Eq. (27) and Eq. (28) ) and k oc ■ Actually, R 2 reflects the 

integration over q t ( b t ) in the first diagrams for the SC ( see Figs. 21 and 22 
). In some sense 



The HERA data on diffractive photoproduction of J/ X I> meson [34] give a 
unique possibility to find R 2 . Indeed, (i) the experimental values for the slope 
( see Fig.23 ) are B e i = 4 GeV~ 2 and B in = 1.66 GeV~ 2 and (ii) the cross 




FIGURE 22. The first order SC ex f d 2 b t fl 2 (s,bt) DIS with the proton in the additive quark 
model for the proton . 



































section for J/\V diffractive production with and without proton dissociation 
are equal. Neglecting the t dependence of the upper vertex in Fig.22, we can 
estimate the value of R 2 . It turns out that R 2 ~ 5 GeV~ 2 or in other word it 
approximately in 2 times smaller than the radius of the proton. Using Eq. (27) 
and this value for R 2 we obtain that A o jo « 40%. 


B Weak SC ( Donnachie - Landshoff approach ): 


Let us ignore everything that I have discussed in the previous section about 
the size of the SC. Let us pretend that we have not learned Eq. (78) - Eq. (81) 
and let us try to estimate the value of the SC using only experimental data. In 
my opinion this was the logic of the Donnachie - Landshoff approach. Indeed, 
they achieved a good description of the experimental data on a tot and a d and, 
therefore, they could expect that the value of SC would be small. However, they 
had to introduce the SC to describe the t- dependence of the elastic cross section. 
Indeed, the experiment shows a clear structure in t dependence - a minimum at 
t « 1.3 — 1.4 GeV 2 ( see Fig.23 ). Everybody knows that such a structure is very 
typical for interaction in optics as well as for interaction of any waves with the 
target with definite size. However, the Pomeron exchange does not reproduces 
this typical diffractive structure. The reason for this is clear since the Pomeron 
corresponds to predominantly inelastic production and strictly speaking the elastic 
cross section should be considered as small in pure Pomeron approach. The t 
structure appears in the Pomeron approach only due to the SC. To understand 
this let us consider a model with one and two Pomeron exchanges. Namely, this 
model corresponds Eq. (78) - Eq. (81). For purpose of simplicity we assume that 
the asymptotic bt - dependence of the Pomeron exchange give by Eq.(69) holds also 
for medium energies. 

Eq. (79) can be rewritten in the form: 


bt) 


ob 



\ 

4 7 tR 2 (s) 



Returning to t representation we have 


A(s,t) = 


(90) 


(91) 


o' o 


JL \ a p (°) e'% 1 1*1 

s o 


167T B el 


Jh)2a P (0) —1 e -%Ut| 

s 0 




where B e i is the t slope of the clastic cross section B e i = h . One can see that at 
sufficiently large —t— the second term in Eq. (91) becomes more important and 
at \t\ = |to| the amplitude has a zero which, actually, manifests itself as a minimum 







FIGURE 24. a e i(s, bt) versus bt in Donnachie - Landshoff approach [35] The values of energies 
are the same as in Fig. 15. 

in the differential cross section. Taking |t 0 | ~ 1.3 GeV 2 one can find A ~ 0.04. 
Donnachie and Landshoff fitted the data using a formula of Eq. (91) - type [35]. 

It should be stressed that this small amount of SC gives about 3 - 5 % contribution 
to the total cross section and, therefore, can be neglected. It is interesting also to 
see how SC helps to preserve unitarity. Fig. 24 shows the ^-dependence of a e i in 
the D-L approach with the SC. One can see that the unitarity will be violated only 
at a/s k 5 TeV (see Fig.24 ). 

C Minimal SC ( Eikonal approach) : 

We assume in the Eikonal approximation, that the hadron states are diagonal 
with the strong interaction. As we have discussed ( see Eq. (54) - Eq. (59) ), in 
this case a SD = 0. Therefore, the Eikonal approach has a well defined precision, 

S P 

namely, we cannot trust the Eikonal approach within the accuracy better than ([—. 
Fig.17 shows that this estimates give errors of the order of 17 - 12 %. This is an 
improvement in comparison with the single Pomeron exchange or with the D-L 
approach with their weak SC, but it is not a perfect approach. 

Let me list here all pluses and minuses of this approach. I hope, that the honest 
discussion of them will lead us to a better understanding of the main theoretical 
problems in high energy scattering. 

9 S l 

1. The Eikonal approach is the first thing that we have to do to comply the 
unitarity constraints ( see Eq. (37) ). Indeed, this approach is the consistent 
way to take into account two main processes: (i) the inelastic interaction with 



Central Production 


0 


Y = ln(S/|l ) 


Diffractive Dissociation 


0 


Y = ln(S/|T) 


FIGURE 25. Structure of the final state in the Eikonal approach. 


the typical Pomeron - like event with large multiplicity of produced particles; 
and (ii) the elastic cross section. The structure of the typical structure of the 
final state is pictured in Fig.25. 

2. The Eikonal approach gives what we expect from the wave picture of the inter¬ 
action. Everybody knows that light ( or any other wave process ) interaction 
with the target gives always a white spot in the middle of the shade. There¬ 
fore, we cannot even dream to be consistent with the quantum mechanique 
without the Eikonal approach. Of course, you can calculate the elastic cross 
section in the D-L approach considering the elastic cross section as a small per¬ 
turbation. However, we have to check how small the calculated elastic cross 
section using the unitarity constraints of Eq. (37) and reconstructing from it 
the total and inelastic cross sections. Fig. 16 shows that the correction are 
large and changes the behaviour of G in (s,b t ). 

3. The Eikonal approach is certainly the simplest way of calculating of SC and 
can be used for the first estimates of how essential they could be. The simplest 
is not always the best and we will show below that the Eikonal approach has 
a serious pathology. 


4. The Eikonal approach describes well the experimental data on total and elastic 
cross sections ( see Ref. [36] and Figs. 11 and 17 ). 



1. In the parton model the Eikonal approach means that only the fastest parton 
















interacts with the target. This is a very unnatural assumption and, in my 
opinion, it is a great defect of this approach. 

SD 

2. In the Eikonal approach the errors are of the order of ((— m 10 — 17% and 
they are not small. 

3. The Eikonal approach cannot explain the CDF data on double parton inter¬ 
action and it give a e ff ~ 30 mb ( see Eq. (88) ) which is in about two times 
large that experimental value ( see Eq. (85) - Eq. (89) ). 

4. The intristic problem of the Eikonal model is the energy behaviour of a SD at 
large masses M 2 oc s. We will show below that in this region of mass the 
diffractive cross section a SD oc s Ap and it exceeds the value of the total cross 
section a to t oc Ins. This is our payment for simplicity and unjustified ( adhoc 
) assumption that only the fastest parton interacts with the target. 

The main assumption of the Eikonal approach is the identification the opacity 
D(s, bt) with the single Pomeron exchange, namely 


D(s, b t ) = POM(s , b t ) ; 


POM(s,b t ) = ; 

R 2 (s ) = i?g + 2 a' P ln(s/s 0 ) ; 

(ip'- = Si(«A = 0) . 


7rR 2 (s) ' SQ 


(92) 

(93) 

(94) 

(95) 


Eq. (92) is certainly correct in the kinematic region where fl is small or, in 
other words, at rather small energies or at high energies but al large values of 
bt . Therefore, the Eikonal approach is the natural generalization according to s- 
channel unitarity of the single Pomeron exchange. Using Eq. (37) and Eq. (92) - 
Eq. (94), one can obtain ( see Ref. [36] for details ): 


c not = 2 7T R 2 (s) { ln(z/(s)/2) + C - Ei(-u/ 2)} ; 
a in = it 7? 2 (<s) {ln(z/(s)) + C - Ei(—v) } ; 

^ei(^) &tot(,s) *^m(^) i 


D I R 2 ( s ) hypergeom(\l ,l,l,l,[2,2,2l,—W2) 

B = -V-lt-0 = -j- ,.„ 0 ”([l,l],[2.2].y/2) 


(96) 

(97) 

(98) 

(99) 


Here, the notations of Maple were used for the generalized hyper geometrical func¬ 
tions. Figs. 11 and 17 show that the Eikonal approach is able to describe the 
experimental data. 

The Eikonal approach gives a very simple procedure how to calculate the survival 
probability of LRG processes ( see our high energy glossary ). Indeed, in this 
approach the general formula of Eq. (68) can be simplified and it has a form: 


< S 2 > = 


/ d\S H (b t ) P(s, b t ) f d 2 b t S H (b t ) e n % 6 -*) 


I d?b t S H (bt 


f d 2 b t S H (bt) 


( 100 ) 



where Snipt) is the profile function for the “hard” processes for two jets production. 

For all numerical estimates we took a Gaussian parameterization for Sn{bt) = 
,2 

r h . In the framework of this simple parameterization all integrals can be 
taken analytically and we have 

2 ^ _ a(s)7[a(s),i/(s)] 

[i/(s)]“W 



with 


B 2 (s) interaction area for “ soft ” collisions 

K 2 h interaction area for “ hard ” collisions 


( 102 ) 


7 [a, u\ is incomplete gamma function. 

Fig.26 shows the Eikonal model prediction for < S 2 > which are in a good 
agreement with the DO data and reproduce the energy behaviour of the LRG rate( 
see Fig.9-b ). 

The key problem in the Eikonal approach is diffraction production which has not 
been included in our unitarization procedure. This is a more complicated way to 
express that DD is equal to zero in the Eikonal approach. However, we can treat the 
diffraction dissociation as a perturbation considering that its cross section is much 
smaller than the clastic one. Fig.17 shows that such an approach cannot not be very 
good since « 50%. The first step is to consider the triple Reggeon diagram 
of Fig.l 8 -b type Appp. The second one is to take into account the rescatterings 
of the fast partons which suppress the value of Appp. Actually, the resulting cross 
section has the same form of the survival probability for LRG processes since the 
DD is the simplest LRG process: 

A/f ^ r 

dM2 U=o = J d 2 b t A P pp(s, M 2 ) P(s,b t ) = (103) 


where 


cr n 




Gppp x 




o(s) 7 [a(s), us 
j/]“(s) 


i 


a(s) = 



(104) 


with = 2 Bf n + r'oj + 4 «p(0) ln(jgy). Bf n denotes the radius of the 

proton - Pomeron vertex while rf x is the radius of the triple Pomeron vertex [37]. 
B 2 (s) = 4 i?oi + 4 a' P In (s). 

Eq. (103) describes experimental data much better than the single Pomeron 
exchange after adding a secondary Reggeon contributions. It shows that we are on 
the correct way. The most important feature of Eikonal approach description is 





the flattering of the energy dependence. However, in Fig.17 we plot the integrated 
cross section over M 2 defining the region of integration 4 < M 2 0.3 Vs 1 GeV 2 . the 
experimental cr SD was integrated in the region 4 < M 2 0.05 s. We cannot trust 
our calculations for M 2 oc s since 


as 


U > i 


2 { 1 


In M 2 
In s 


(105) 


Using the asymptotics for the incomplete gamma function we obtain that 


M 2 da SD 
dM 2 


( M 2 ) Ap . 


(106) 


Integrating Eq. (14) over M 2 up to M 2 = ks one obtain 

rks rj A/f 2 1 

V SD OC J -jp(M 2 ) Ap = —( ks) Ap > a tot oc Ins. (107) 


Eq. (16) shows that there is no way to incorporate the diffractive production in the 
Eikonal model since the “small” parameter of such an approach cannot be small. 

This example is the most transparent illustration of the difficult problem in taking 
into account the SC - the description of high mass diffraction. 


D Correct SC ( ? ) 

Certainly, we do not know an answer to the question what is a correct procedure 
to take into account SC. Here, I discuss two models for SC which both are better 
than Eikonal approach. 

Quasi-Eikonal Model. 

The idea of this model is very simple: to take into account the diffractive dissoci¬ 
ation in the same Eikonal way as clastic rescatterings. Doing so, we assume that we 
have two wave functions which are diagonal with respect to the strong interactions: 
dq and d/ 2 . The general Eq. (55) has the form 


hadron ~ Pl X ^ 1 + P2 d / 2 j (108) 

with condition: + /3| = 1, which follows from the normalization of the wave 

function. The wave function of the diffractively produced bunch of hadrons should 
be orthogonal to d J hadron and looks as follows: 

d >d = — /^du + Pi d/ 2 • (109) 

The elastic and single diffraction amplitude in this model can be rewritten, using 
Eq. (56), Eq. (108) and Eq. (109) in the form: 

a e l(s,b t ) = < d 'final T hadron > = P\ M + P\ A 2 ; (110) 

aSD(s,b t ) = < d 'final d 'd > = Pi P 2 { ^2 — ^1 } • (HI) 



FIGURE 27. The Pomeron exchanges in the Quasi-Eikonal model. 


For amplitude A n (s,b t ) we have a unitarity constraint of Eq. (60) and therefore: 

A n {s,b t ) = i{ 1 - e } ; (112) 

G n (s,bt) = 1 - . (113) 


Assuming the simplest eikonal form for Q n ( see Eq. (92) - Eq. (95) ), namely 


Hn(s, bt) 


&0 n 

nR 2 n (s) 




5 


we have the final answer: 


otot = 2 f d?b t {0j [ 1 — e 2 ] + (1 - ft ) [ 1 - e s 2 ] } ; 
a a = J d 2 b, [ft [1 - e-%- ] + (1 - ft ) [ 1 - e} ] 2 ; 
a SD = /<f 2 6,ft(l - ft)[e-T? - e-T] 2 . 


(114) 


(115) 

(116) 
(117) 


It is obvious from Eq. (115) - Eq. (117) that we introduce three new parameters to 
describe SC including the processes of diffractive dissociation: /?| cr 2 o and /?i. They 
are correlated with the total cross section of the diffractive dissociation and its in¬ 
dependence. At first sight it is strange that we have to introduce three parameters 
not two. The third parameter cr 0 2 is closely related to the cross section of “elastic” 
rescattering of the hadrons produced in the diffractive dissociation. The meaning 
of all parameters is clear from Fig. 27. 

We can rewrite Eq. (115) - Eq. (117) in more convenient form, using the experi¬ 
mental fact that the exponential parameterization works quite well for differential 
cross section of DD. Indeed, we can assume that 


0 2 - Hi 


11 d (Sj bt) 


o'DO / s 

7T R 2 d(s) So 


Ap e & D (s) 


(118) 


while for Oi we use the parameterization of Eq. (114). In such a parameterization 
we have a very simple form for a e i ands an, namely 


o e i(s,6 ( ) = 1 - e a 1 + (l-ft)e 2 1 { 1 - e ? } ; 
a D (s,b ,) = ft (1 - ft) { 1 - e~~f~ } e-S 1 . 


(119) 

(120) 



Fig. 28 shows the result of fitting of energy behaviour of different 
ables. The parameters, that were used, are 
1. A P = 0.1 ; 


‘soft” observ¬ 


ers. 


= 0.06 


VD 

W. 


50 

7T Ki 


3 . 02 = 0.5 


4. R\ = 26 + 4ap(0) ln(s/s 0 ) with o;p(0) = 0.25 GeV~ 2 . 

5. R 2 D = 13 + 4ap(0) ln(s/so) with ttp(0) = 0.25 GeV~ 2 . 

The only parameter which we have to comment is a sufficiently big difference in 
for elastic and diffractive scattering. We think, that it happened because of the 
fact that we took effectively into account the integration over produced mass in our 
parameterization. 


“Fan” diagrams (Schwimmer resummation ). 

The technical problem that we are going to solve here is the resummation of 
all “fan” diagrams (see for example Fig. 18-b and Fig. 29-a). it is interesting 
to mention that there is a problem which can be solved by such a resummation, 
namely, the interaction of a fast hadron with heavy nuclei ( see Refs. [38] [39] ). 
Indeed, the vertex for Pomeron - nucleus interaction is proportional to g n A, where 
gisr is the Pomeron-nucleon vertex and A is the atomic number. Therefore, the 
“fan” diagrams give you the maximal contribution ( oc (gvGApA)" where G 3 p is 
the vertex for triple Pomeron interaction ). Neglecting all momentum transferred 
dependence in nucleon - nucleon scattering since R% = R$ + ap(0) ln(s/so) -C 
R\ we can write a simple equation which sums all “fan” diagrams. The graphic 
form of this equation is presented in Fig.29-b and it has the following analytical 
form for p + A interaction at Y = ln(s/so): 


G P (Y,b t ) = g 2 N AF A (b t )e ApY - g N G 3P C dy e Y G 2 P (y,b t ) . (121) 

Jo 

Taking derivative with respect to Y and using Eq. (121) we can reduce this 
equation to simple differential equation: 


dG P (Y.b t ) 

dY 


A P Gp(y.b t ) ~ 7 G 2 P (Y,b t ) , 


( 122 ) 


where gamma = g n G 3 p . The solution should satisfy the following initial condition 
which is obvious directly from Eq. (121): 


Gp(Y,b t )\ Y =o — 9 2 NAF A (b t ) , 


(123) 


where F A (b t ) directly related to form factor of a nucleus F(t) = / d 2 b t e l ^ bt F A {b t ) 
with q 2 = —t. 

One can obtain the solution: 


G P (Y,b t ) = 


g 2 N AF A {b t ) e 


A P Y 


1 + g 2 N AF A (b t )£ { e ApY - 1} ' 
One can see that this solution has a nice properties: 


( 124 ) 







( 125 ) 


1. If AF A (b t )e ApY > 1 G P 

2. If AF A (b t )e ApY « 1 G P 

3. the equation 

g 2 N AF A (b 0 ) -j-e A ' Y = 1 

lAp 

defines the value of bt = bo, which is the border between the above two con¬ 
ditions. Namely, for b t < b 0 the first one holds, while for b t > bo the form 
factor is so small that we are in the region of the second constraint. 

4. Talking for simplicity the exponential form for form factor F A (bt) = 

i2 

_ 

—kre r a we obtain the solution to Eq. (125): 

’ ku a 

bl = R\ A P Y ln(g 2 N AF A (b 0 ) )^~) . (126) 

Ap 

5. From Eq. (126) we can obtain the asymptotic formula for the total cross sec¬ 
tion: 

cr tot (nA) =>• 27t6q oc A$ A + pY ln(AY) (127) 

6. The cross section of the diffractive dissociation turns out to be constant at large 
M 2 ~ s. Therefore, this approach has no such difficulties as Eikonal and/or 
Quasi-Eikonal Model. The reason for this is quite simple: In the Eikonal 
model we took into account only the rescatterings of the fastest parton (see 
Fig.l8-a ) in the parton cascade. The “fan” diagrams describes the interaction 
of all partons with the target as it can be seen from Fig.l8-b. 

7. Eq. (127) gives a nice illustration that we can have an increase in the radius 
of interaction without assuming any slope for Pomeron trajectory. 

The “fan” diagrams gives a rather selfconsistent approach which heal a lot of dif¬ 
ficulties but the main problem is to justify this approach. I gave one example - 
hadron - nucleus collision- for which the problem can be reduced to summation of 
the “fan” diagrams. You will see more in this lecture. 

To finish with hadron - nucleus interaction we have to treat Gp as opacity in the 
eikonal formula Q — Gp (see Eq. (92) - Eq. (95) ). 

II “H A R D” POMERON 


Ap 
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VI 


EVOLUTION EQUATIONS AT X —> 0 
Parton cascade and evolution equations: 


A 


As it has been mentioned ( see Eq. (15) - Eq. (14) ) “hard” Pomeron is nothing 
more than the high energy behaviour of the scattering amplitude at short distances 
that follow from the pQCD. In spite of the fact that solution of the evolution 
equations gives nothing like the Regge pole - Pomeron, the parton cascade picture 
( see Fig.l ) is still correct for interpretation. Let me derive for you Eq. (14) and 
show you that that simple equations (Eq. (1) - Eq. (5) ) work quite well for the 
“hard” Pomeron too. 

The probability to emit one extra gluon in QCD is equal to 


Pi 


N c as(ku) dxj dk 2 t 

IT Xi kl 


(128) 


where x^ is a fraction of energy carried by gluon “i” and k t i is its transverse mo¬ 
mentum. Factor is nothing more than the invariant volume for emission 

( d 3 ki/Ei ) and the extra k 2 comes from the gluon propagator. You can check 
that Eq. (128) gives the only possible dimensionless combination. Recall, QCD is 
dimensionless theory. 

To calculate the probability to emit n ( W n ) gluons we need to find the product 
of probabilities, namely 


w n = n Pi, (129) 

i =o 

and to find the probability to emit arbitrary number of gluons we have to sum 
Eq. (130) over n: 

oo oo n 

W = £ W n = £ n p ‘ (!30) 

n= 0 n=0 7=0 

^ A A N c a s (kl) dxj dk | 

2^ II p ■ 

71=0 7=0 n h ti 

To obtain the deep inelastic structure function we have to integrate over all possible 
emission with the conditions: Xi > Xb and k\ < Q 2 where Xb = — and q 2 
is the virtuality of the probe. The largest contributions we can obtain from the 
phase space in which both integrals with respect xi and k 3 ti give logs ( Iti(1/xb) 
and InQ 2 , respectively). This phase space is the following: 


» ii » i 2 > ••• » » x i+ i S> ... 3> x n » xb ; 

Q 2 » k 2 n > kl n _ x > ... > k 2 n > Ql . 


(131) 

(132) 



Therefore, we have a strong ordering in transverse momenta and fractions of 
energy. Taking into account Eq. (131) and Eq. (132), we can easily rewrite ( in the 
simplest case, neglecting k 2 dependence of QCD coupling constant ) 



N c a s (k 2 i) dxj dk 2 t 

Tl Xi k 2 ti 


_ 1 _ { N c a s m HQ 2 IQl) h1/xb) y (133) 


Therefore, the structure function can be found iust summing all contributions of 
Eq. (133): 


xG(x, Q 2 ) 


OO 


E 

n=0 


L n 

n\ 2 


(134) 


where scale L is equal to 

L = N(: ° s ln{Q 2 /Ql ) ln(l/ xb) = £ y (135) 

7r 

where £ = Vcas J Cj » ] \n(Q 2 /Q'fy and y = ln(l /xb)- For large Q 2 and small xb 
parameter L 1 and, therefore, only sufficiently large n will be dominant in 
Eq. (135). Since 


2n! = T(2n + 1) 


2 2n j 

—=T(n+ 1) T(n + -) 



we can rewrite Eq. (134) in the form 


xG(x, Q 2 ) 


OO 


E 

n —0 


(2vT) 2n 

2ffi 


e 2VZ 


(136) 


(137) 


On the other hand, using Stirling formula, we can evaluate 

2n! -> (2 n) 2n e~ 2n , 


which, being substituted in Eq. (137) gives the estimate for the average n, namely, 

< n > = 2 VL . (138) 

One can see that Eq. (137) can be written as 

xG(x, Q 2 )|l> i ==>■ e <n> . (139) 


For running coupling constant the only difference is the redefinition of the variable 

« =iii S(S iiiEq ' (135) ' 
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FIGURE 30. Parton cascade and evolution equation ( Feyman diagrams). 


B Why do we use evolution equations? 

The first question is what is the meaning of < n >. To answer this question we 
need to look in the picture of our cascade more careful. Let us assume that the 
number of “wee” partons N is not extremely large or in other words there is only 
one “wee” parton with definite fraction of energy (x{ ) and transverse momentum 
(kti ). In this case only one “wee” parton interacts with the target and all others 
gather together without loosing their coherency. In the Feynman diagrams they 
contribute to the renormalization of mass and coupling constant in QCD. This is 
shown in Fig.30. 

In Fig. 30 one can see that < n > means the average number of cells in our 
“ladder” diagram or, in other words, in our evolution equations. On the other had, 
< n > is the average number of produced jets ( mini - jets ) in the whole region 
of rapidity accessible in the experiment. I found very instructive to plot what we 
know from HERA experiment about the value of < n > ( see Fig.31 ) 

In Fig. 32 < n > is plotted down to x — 1CT 6 and looking at this plot we can 
see, that 

• < n > increases mostly due to cascading rather in Q 2 than in ln(l/x) ; 

• The maximum value of < n > is about 5 ; 

• < n > is the total multiplicity of produced jets and/or minijets. Therefore, the 
average density of produced jets ( minijets) is less than d < n> jdy ~ 0.5. 
























































FIGURE 31. < n > versus x at different values of Q 2 from HERA experiment. 








This number suggests that we need to calculate only the Feynman diagrams 
of up to a 2 s or/and order ; 

• High energy resummation is only a way to get N without calculating < n > . 

The main conclusion from this discussion is very simple: we have to use the 
evolution equations only because it is the only known way how to calculate N 
without evaluation of < n >. However, if we will find a way to calculate < n > 
we will be able to predict our structure functions with better accuracy restricting 
yourselves by calculation only limited number of diagrams. 


VII SC FOR “HARD” POMERON 

The main outcome of our previous discussion is the fact that N 1 at high 
energies. It happens both for “hard” ( see Eq. (14) and/or Eq. (138)) and “soft” 
processes. Therefore, in the collisions at high energy the new interesting system of 
parton has been created with large density of partons. For the “hard” processes 
all these partons are at the short distances where the coupling QCD constant is 
small enough to be use as a small parameter. However, we cannot apply for such 
a system the usual methods of perturbative QCD because the density of partons 
is large. In essence, the theoretical problems here are non perturbative, but the 
origin of the non perturbative effects does not lie in long distances and large values 
of as, which are typical for the confinement region. 


A Which parton density is large? 


The qualitative estimates of which density is high, we can obtain from the s- 
channel unitarity [11] which we can use in two different form: 

• Vtoti'fp) — a em2vri? 2 , where R is the size of the target, a em is the fine 
structure constant ( see Ref. [43] for explanation why we have a em ). This 
constraint can be rewritten as [11] [12] [40] 


K 


3 7T asxG(x, Q 2 ) 
Q 2 R 2 


(140) 


• a DD ('y*P ) < Vtot(j*p ), where DD stands for diffractive dissociation. This 
inequality leads also to Eq. (140). 

Therefore, 

1. If k is very small (k € 1 ), we have a low density QCD in which the parton 
cascade can be perfectly described by the DGLAP evolution equations [5]; 




ln(Q 2 /GeV 2 ) 

FIGURE 32. Contours for < 7 > and k for the GRV’9f gluon density and HERA kinematic 
region. 


2. If k < 1, we are in the transition region between low and high density QCD 
. In this region we can still use pQCD, but have to take into account the 
interaction between partons inside the partons cascade; 

3. If k > 1, we reach the region of high parton density QCD, where we need to 
use quite different methods from pQCD 

First, we want to make a remark on parton densities in a nucleus. Taking into 
account that for nucleus Ra = Rn x As and xGa(x, Q 2 ) = AxGn(x,Q 2 ), we 
can rewrite Eq. (140) as 


k a 


3 7 r as A xGn(x, Q 2 ) 

Q 2 R% A% 


A 3 kn A 1 • 


(141) 


Therefore, for the case of an interaction with nucleii, we can reach a hdQCD region 
at smaller parton density than in a nucleon . 

Fig.32 gives the kinematic plot (x, Q 2 ) with the line k — 1, which shows that the 
hdQCD effect should be seen at HERA. 


B Two different theoretical approaches 

To understand these two approaches we have to look back at Fig.l which gives 
you a picture of high energy interaction in the parton approach. In Eq. (1) N 
is the flux of partons. The main idea of SC is the fact that this flux should be 
renormalized in the case of high parton density. Indeed, if the number of “wee” 






partons is so large that, let say, two partons have the same kinematic variables : Xi 
and k t i, we do not need to take into account the interaction of these two partons two 
times. Let me recall, that the total cross section is the number of interaction but 
not the number of partons. It means that if the first of two partons has interacted 
with the target the total cross section does not depend upon the fact did the second 
interact or not. It is obvious that the probability to have two partons which could 
interact with the target semalteneously is equal to 


2 nR 2 ( Xi ) ’ 

where irR 2 is the area which is populated by partons with the fraction of energy 
Xi. Therefore, we have to renormalize the flux 



RENORMALIZED FLUX = N, 


REN 


= N x {1 - P 2 } . 


(143) 


One can see that this renormalized flux gives the total cross section in the form: 


^ _ ^hp 

&tot. — otot 
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HP \2 
tot ) 


vr K 2 


(144) 


which is a Glauber formula for SC in the limit of small SC. 


However, when N 1, we have to 

Tr , r , ,, , ,. change our approach completely from the 

If IV « 1, we expect that the , , i , 

parton cascade to one based on 

semiclassical field approach, since due to 

the uncertainty principle ANAcj) « 1, 

we can consider the phase as a small 

parameter. Therefore, in this kinematic 

region our magic words are: 

• Semi-classical gluon fields; 

• Wiezsacker-Williams approximation; 

• Effective Lagrangian for hclQCD; 

• Renormalization Wilson group 
Approach. 

It is clear, that for N & 1 the most natural way is to approach the hclQCD 
looking for corrections to the perturbative parton cascade. In this approach the 
pQCD evolution has been naturally included, and it aims to describe the transition 
region. The key problem is to penetrate into the hclQCD region where k is large. 
Let us call this approach “pQCD motivated approach ”. 

For N 3> 1, the most natural way of doing is to use the effective Lagrangian 
approach, and remarkable progress has been achieved both in writing of the explicit 
form of this effective Lagrangian, and in understanding physics behind it [13]. The 
key problem for this approach was to find a correspondence with pQCD. This 
problem has been solved [45]. 

Fig.33 shows the current situation on the frontier line in the offensive on hclQCD. 


renormalization of the flux 
will be small, and we use an 
approach with the following 
typical ingredients: 

• Parton Approach; 

• Shadowing Corrections; 

• Glauber Approach; 

• Reggeon-like Technique; 

• AGK cutting rules. 
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FIGURE 33. The current situation in the battle for hdQCD. 


C The picture of interaction 

To understand the picture of interaction in the region of small x it is better to by 
examine the parton distribution in the transverse plane (see Fig.those partons with 
size r p ~ A. At x approx 1 a few parton are distributed in the hadron disc. If we 
choose Q such that r 2 -C K\ then the distance between partons in the transverse 
plane is much larger than their size, and we can neglect the interaction between 
partons. The only essential process is the emission of partons,which has been taken 
into account in QCD evolution. As x decreases for fixed Q 2 , the number of partons 
increases, and at value of x = x cr , partons start to populate the whole hadron disc 
densely. For x < x cr the partons overlap spatially and begin to interact throughout 
the disc. For such small x values, the processes of recombination and annihilation 
of partons should be as essential as their emission. However, neither process is 
incorporated into any evolution equation. What happens in the kinematic region 
x < x cr is anybody’s guess. We suggested that parton density saturates, i.e. the 
parton density is constant in this domain. 


D The GLR equation 

The first attempt to take into account the parton - parton interaction in the 
pQCD motivated approach was done long ago [11]. It was based on the simple idea 
that there are two processes in a parton cascade (see Fig.l) : (i) the probability 





of the emission of an extra gluon is proportional to as p where p is the density of 
gluon in the transverse plane, namely 


P 


xG(x , Q 2 ) 


71 R 2 


(145) 


and (ii) the annihilation of a gluon, or in other words a process in which the 
probability is proportional to p 2 x cr annihi i ation . <T annihi i ation can be estimated as 
aannihilation oc as r r 2 , where r is the size of the parton produced in the annihilation 
process. For deep inelastic scattering, r 2 oc Therefore, in the parton cascade 
we have 


Emission (1 —> 2 ) : probability = p emisslon oc as p ; (146) 

Annihilation (2 —> 1 ) : probability = panmhiiation ^ pi 2 . (147) 

At x ~ 1 only emission of new partons is essential, because p •Cl and this 
emission is described by the DGLAP evolution equations. However, at x —> 0 the 
value of p becomes so large that the annihilation of partons becomes important, 
and so the value of p is diminished. The competition of these two processes we 
can write as an equation for the number of partons in a phase space cell ( Ay = 
A ln(l/x) A \n(Q 2 /Ql) ): 
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or in terms of the gluon structure function xG(x, Q 2 ) 
d 2 xG(x, Q 2 ) as N ( 
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(149) 


This is the GLR equation which gave the first theoretical basis for the consideration 
of hdQCD. This equation describes the transition region at very large values of Q 2 , 
but a glance at Fig.32 shows that we need a tool to penetrate the kinematic region 
of moderate and even small Q 2 . 


E Glauber - Mueller Approach 

We found that it is very instructive to start with the Glauber approach to SC. 
The idea of how to write the Glauber formula in QCD was originally formulated 
in two papers Ref. [41] and Ref. [40]. However, the key paper for our problem 
is the second paper of A. Mueller, who considered the Glauber approach for the 
gluon structure function. The key observation is that the fraction of energy, and 
the transverse coordinates of the fast partons can be considered as frozen, during 
the high energy interaction with the target [41] [40]. Therefore, the cross section 



of the absorption of gluon((S'*) with virtuality Q 2 and Bjorken x can be written in 
the form: 
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where z is the fraction of energy which is carried by the gluon, T'j is the wave 
function of the transverse polarized gluon, and cr N (x,r 2 ) is the cross section of the 
interaction of the GG- pair with transverse separation r t with the nucleus. Mueller 
showed that Eq. (10) can be reduced to 
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and prohle function S(b t ) chosen in the Gaussian form: S(b 2 ) = . 

Obviously, the Mueller formula has a defect, namely, only the fastest partons 
(GG pairs) interact with the target. This assumption is an artifact of the Glauber 
approach, which looks strange in the parton picture of the interaction. Indeed, in 
the parton model we expect that all partons not only the fastest ones should interact 
with the target. At first sight we can solve this problem by iteration Eq. (151)( see 
Ref. [42] ). It means that the first iteration will take into account that not only the 
fastest parton, but the next one will interact with the target, and so on. 


F Why equation? 


We would like to suggest a new approach based on the evolution equation to 
sum all SC. However, we want to argue first, why an equation is better than any 
iteration procedure. To illustrate this point of view, let us differentiate the Mueller 
formula with respect to y — ln(l/x) and £ = InQ 2 . It is easy to see that this 
derivative is equal to 


d 2 xG(x, Q 2 ) 
dy <9£ 



e 






(152) 


The advantages of Eq. (152) are 

• Everything enters at small distances ; 

• Everything is under theoretical control ; 



• Everything that is not known ( mostly non perturbative ) is hidden in the 
initial and boundary condition. 

Of course, we cannot get rid of our problems changing the procedure of solution. 
Indeed, the non-perturbative effects coming from the large distances are still impor¬ 
tant, but they are absorbed in the boundary and initial conditions of the equation. 
Therefore, an equation is a good ( correct ) way to separate what we know ( short 
distance contribution) from what we don’t ( large distance contribution). 


G Equation 


We suggest the following way to take into account the interaction of all partons 
in a parton cascade with the target. Let us differentiate the ^-integrated Mueller 
formula of Eq. (151) in y = ln(l/x) and £ = ln(Q 1 2 3 /Qo)- ft gives 

d 2 xG(y, £) _ 2 Q 2 K 2 f n , w ( , ^ 

—QyQc — — — ^2 — { C tQ )) + Ei(kg{x ,Q ))| , (153) 

where kg(x, Q 2 ) is given by 

k(x,Q 2 ) = 2Q2R2 xG ( x ’ ^ 2 ) • ( 154 ) 

Therefore, we consider kq on the l.h.s. of Eq. (153) as the observable which is 
written through the solution of Eq. (153). 

Eq. (153) can be rewritten in the form ( for fixed as ) 


d~K G {y,£) dn G {y,£) 

dyd^ dy 

{ C + ln(/c G ) + E 1 (k g )} = F(kg) 

7T 

This is the equation which we propose. 

This equation has the following desirable properties: 


(155) 


1. It sums all contributions of the order (asyC ,)", absorbing them in xG(y,£), 
as well as all contributions of the order of n n . Therefore, this equation solves 
the old problem, formulated in Ref. [11], and for N c —> oo Eq. (155) gives the 
complete solution to our problem, summing all SC ; 

2. The solution of this equation matches with the solution of the DGLAP evolu¬ 
tion equation in the DLA of perturbative QCD at k —» 0 ; 

3. At small values of k ( k < 1 ) Eq. (155) gives the GLR equation. Indeed, 
for small k we can expand the r.h.s of Eq. (155) keeping only the second 
term. Rewriting the equation through the gluon structure function we obtain 
the GLR equation [11] with the coefficient in front of the second term as 
calculated by Mueller and Qiu [12] ; 
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FIGURE 34. The estimates of SC based on Eq. (155). 

4. The first iteration of this equation gives the Mueller formula ( see Ref. [40]) ; 

5. In general, everything that we know about SC is included in Eq. (155) ; 

H The size of SC 

In Fig.34 we plot our calculation for the DGLAP evolution equation, for the 
Glauber - Mueller approach, for the GLR equation and for the new evolution equa¬ 
tion. We can conclude, that: 

1. SC even in the Glauber - Mueller approach are essential for the gluon structure 
function ; 

2. The Glauber - Mueller approach considerably underestimates the value of SC ; 

3. The GLR equation leads to stronger SC than the solution to Eq. (155) ; 

4. The new evolution equation does not reproduce the saturation of the gluon 
density in the region of small x which the GLR equation leads to . 

We firmly believe that the new evolution equation gives the correct way of eval¬ 
uation of the value of SC. However, the difficult question arises: why the SC have 
not been seen at HERA? Our answer is: 

1. The value of SC is rather large but only in the gluon structure function while 
their contribution to F 2 is rather small [42] (see Fig.34 ), and cannot be seen on 
the background of the experimental errors ; 

2. The theoretical determination of the gluon structure function is not very 
precise and we evaluate the errors as 50% [44] and the SC in xG is hidden in such 
large errors ; 

3. The statement that SC have not been seen is also not quite correct. Our 
estimates [46] show that the contribution of SC is rather large in the F 2 slope 





( - t’,, -z’j) 



FIGURE 35. The space - time structure of the parton cascade. 
and incorporating the SC one can describe the recent experimental data [47] on 

dHQVcfi) (see Ref - t 40 ] for detail )■ 

I Effective Lagrangian Approach. 

Deeply in my heart, I firmly believe that there will be a bright future for the 
effective Lagrangian approach in which one combines the physics of hdQCD, with 
the formal methods of the quantum held theory, and gives a way to incorporate the 
lattice calculation for the non-perturbative observables. However, the parameter or 
better to say a new scale Ql(x) which is the solution of the equation k(x, Qq(x)) = 
1 only can be found in the pQCD motivated approach, since the effective Lagrangian 
was only derived and justified assuming this new scale. 

Here, 1 would like to show what kind of physics have been used to obtain the 
effective Lagrangian. The space - time picture of the parton cascade is given in 
Fig.35. The main and very important feature of this cascade is that in the region 
of low x ( in so called leading log (1/x) approximation of pQCD ) is the fact that 
a parton with higher energy in the cascade lives much longer than a parton with 
smaller energy. Let us look for the interaction of the parton emitted at time t 2 . 
This parton lives much shorter time than all partons emitted before. Therefore, all 
these emitted before partons will have enough time form a current which depends 
on their density. We assume that this density is large and because of this fact 
our current actually is a normal classical current. Finally, the Lagrangian of the 
interaction of the parton emitted at time t' 2 can be written as 

L(p) + jp(p) • + L{A) , 

where A is the held of a parton emitted at t 2 . However, we can consider a parton 
emitted at t — t 3 and include the previous one in the system with density p. The 

























form of Lagrangian should be the same. This is a very strong condition on the form 
of Lagrangian, so called Wilson renormalization group approach. In Ref. [13], this 
approach has been used to obtain the effective Lagrangian. 



Ill THE BFKL POMERON 


VIII THE BFKL PARTON CASCADE AND ITS 

EVOLUTION. 

The BFKL Pomeron [7] is an asymptotic of the scattering amplitude in pertur¬ 
bative QCD in the kinematic region where the log scale L = as \n{l/x) 3> 1 but 
the values of virtualities of incoming particles are moreless the same Qq ~ Q 2 . 
Practically, this means that we calculate the same diagrams as before but in the 
kinematic region where Eq. (138) holds while there is no special requirements for 
transverse momentum integration ( there is no Eq. (139) ). Of course, the BFKL 
Pomeron should match the “hard” Pomeron in the situation when Q 2 gg Q q. 

First thing, that 1 would like to stress is the fact that the origin of the ln( 1/x) 
contributions is very simple. As it has been mentioned the factor d,Xi/xi in Eq. (128) 
comes from the phase space. Therefore, if we write the optical theorem for total 
cross section of two particles with virtualities Q 2 and Qq in the form 

°° f fj f ■ 

Vtot(x,Q' 2 ,Ql) = E Y[—d 2 k u \M 2 ^n(xi,ki)\ 2 (156) 

n =2 J i X i 

on can see that to obtain log(l/x) contribution we can safely put all Xi = 0 in 
amplitudes M n . It gives 


Vtot(x,Q 2 ,Q 2 0 ) 


E /II— x (d 2 k ti \M 2 ^ 2 +n {Qik )\ 2 = 
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Deriving Eq. (157) we used the ordering in Xi of Eq. (138). 

Assuming that the inte¬ 

gral over k t i in Eq. (157) gives / d 2 hu \M 2 ^ 2+n (0, ki)\ 2 = we can easily 
obtain that 


om( x , Q 2 ,Ql) 
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where Ul oc as- The result of Eq. (158) can be rewritten in the form a to t oc e <n> 
with < n >= ul ln(l/x). Note, that factor in front can be obtained just from 
dimension of the total cross section and the fact that the scattering amplitude 
in the region of Q 2 ~ Ql should be symmetric with respect to Q 2 —> Q q and 

Qg -»• Q 2 - 



It is interesting to point out that the simple feature of QCD,namely, the fact that 
we have a dimensionless coupling constant in QCD. In such theory each emission 
leads to the change of the In k t of a parton by some value ( d ) : 

Ain (kl/Ql)=d. (159) 


After n- emissions the total change of the transverse momentum will be equal to; 


< In (k ti /Q q) > I after n 


emissions 
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(160) 


It means that in average for the scattering amplitude we have 

< In 2 (Q 2 /Ql) > = d < n > — dtUL ln(l/a;) = 4 D ln(l/x) . (161) 


Therefore, the transverse momentum dependence of the scattering amplitude is a 
distribution with < ln(Q 2 /Ql ) > given by Eq. (160). 

Collecting everything that has been discussed we have 
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Eq. (162) gives the answer. The values of three constants: <j 0 = ck| , wi oc 
as and D oc 0:5 have been found in the BFKL approach to perturbative QCD 
[7]. It should be stressed that we have used heavily a fact that out theory has 
no dimension scale. Actually, it has because the QCD coupling depends on the 
transverse momentum. Therefore, strictly speaking, Eq. (162) is valid only for 
fixed QCD coupling. 


IX THE BFKL EQUATION. 

A cks(Qo) fixed. 

This section will be a bit more formal than everything that was before. Never¬ 
theless, I would like to share with you a beautiful result, that the “soft” Pomeron 
might be the BFKL Pomeron with running QCD coupling + confinement. 

The natural function for which the BFKL equation can be written is the gluon 
density (j)(y , q 2 ) which is closely related to the gluon structure function. 

xG(x,Q 2 ) = J ois(q 2 ) 4>(y, q 2 ) , (163) 

where we introduce from the beginning a new symmetric variable y = ln^^. In 
the leading order y is equal to y = ln(l/x) but in the next to leading order this 
variable turns out yo be much more convenient. 

The BFKL equation [7] has a form: 



CO 


G Jq- - q) 


\ 



FIGURE 36. The BFKL equation. 


d(p(y,q 2 ) 

dy 


(164) 


5 {2 \(q - q 0 ) + J d 2 q K(q,q)(j)(y, q,q 0 )) - 2 a G (q 2 ) 0 ( y, q, q 0 ) ) , 

where all notations are clear from Fig.36. 

The BFKL kernel has a remarkable property that function 

0 7 ( q2 ) = ( 165 ) 

is an eigenfunction, namely 

J d 2 q' K{q,q')(j) u {q') = u(f) 0 /(g 2 ) , (166) 

with 

w(/) = {*(/) + *(!-/)_ 24/(1)} = (167) 

+ D{f-\) 2 + 0((/-|) 3 ) , (168) 

where 4/(z) = z (r(z) is the Euler gamma-function. For the Mcllin transform 
of 7T 


0 ( 2 /, <? 2 ) 


f*a+ioo 


da; 

2 7rf 


e^G^q 2 ) 


(169) 


the BFKL equation reduces to the simple form: 

















u>GM) = 1 + u(f)GM) , 


( 170 ) 


where 

ra+ioc sj f 

GM 2 ) = / e-« r - r « ) GM) , (171) 

■J a—ioo Z 7X X 

with r — r 0 = ln(Q 2 / Qq) . 

Substituting solution to Eq. (170) G u (f ) = into Eq. (169) and Eq. (171), 

we obtain after integration over u since we can close contour on the pole in u 

ra+ioo rl f 

<p(y,q 2 )= GM) , (172) 

J a—ioo ZTXl 

where G u (f) should be found from initial nonperturbative parton distributions. 
Using the expansion of ou(f) at / —> one can see that we are able to take the 
integral by saddle point method and obtain the answer of Eq. (162) with 

u L = Ncas J Q o ] 4 In 2 ; (173) 

D = Nca ^ Q & 14 C (3) . (174) 


B Running QCD coupling. 

The previous subsection was written mostly because I would like to point out 
that the running QCD coupling leads to a result which is very important since it 
gives the only theoretical way that I know to obtain the Pomeron - a Reggeon with 
intercept bigger than 1. Indeed, let us consider the BFKL kernel to be proportional 
to the running as(q 2 ) = where r 0 = InQl/A 2 ) and r = InQ 2 /A 2 ) . A 2 

is the position of Landau pole in as where QCD coupling constant becomes very 
large. The generalized BFKL equation [48] looks as follows: 


#(jU q 2 

dy 


(175) 


5 (2) ((q - q 0 ) + — f d 2 q' K(q,q , )(p(y, tf,q 0 )) . 
r J 


This equation can be easily rewritten as the differential equation in double Mellin 
transform with respect to y and r ( see Eq. (165) ). 


< ^fi , (^)/') / r\ ( r\ | —fr n 

- = roUconfU) g{u, J) + r o e 70 . 


(176) 


The solution of homogeneous equation ( Eq. (176) without the last term ) can be 
easily found and it has the form ( see Refs. [48] for details ): 



(177) 


9(w,f) 


g(u)e 


- 7? ff Q “confU'W 


Function g(ui) should be specified from initial or boundary conditions. The value 
of /o can be arbitrary since its redefinition is included in function g(u). Unless it 
is specially stipulated f 0 = 0. 

Let us find a Green function of the BFKL equation with running QCD coupling 
( G r (y,r ) ) which satisfies the following boundary condition: 


Gr(y,r ) : 

G r (y,r = r 0 ) = 8(y - y 0 ) 


(178) 


This Green function allows us to find us the solution of the BFKL equation for any 
boundary input distribution G in (y,q 2 = q^) at q 2 = q$ ( r = r 0 ). Indeed, such a 
solution is equal to 

G(y,r) = j dy 0 G r (y,r)G in (y 0 ,q 2 = qjj) . (179) 

Such a Green function is very useful for study of the boundary condition for 
the DGLAP evolution. Using G r (y,r) and Eq. (179), we can investigate the in¬ 
dependence at q 2 = q 2 y We can distinguish two cases with different solutions: 

1. The integral over y 0 depends mostly on properties of input function G m ] 

2. The integral over y 0 is sensitive to the Green function. In this case we can 
claim that the energy behaviour of our boundary condition is defined by the 
BFKL dynamics. 

Therefore, this Green function ( G r (y,r) ) can provide us with an educated guess 
for the energy dependence of the boundary condition in the DGLAP evolution 
equations [5]. 

Substituting in Eq. (177) the expansion of c o con f at / —> ^ ( see Eq. (168) 
) we obtain that integral over / in Eq. (177) gives the Airy function 
^ 5 t r — )• Therefore to satisfy the boundary condition of Eq. (178) 

we have to choose a function g(u>) = AG 1 ((^^)^ [p) — yf r o])- 

Finally [48], G r (y — y 0 , r, r 0 ) is equal to 

G r (y ~ yo,r,r 0 ) = (180) 

nr p+ioo dy^ (v _ vo) 

y ro Ja-ioo 2m Ai((^)5[r° - ^r 0 ]) 

Airy function is an analytical function which has zeros and a singular behaviour 
at large values of the argument Ai(z) | 2> o ; p|>i —> -m-e~ ■ Therefore, we have 

two possibility to take the integral over lu (see Fig.37, which shows the structure of 
singularities in uj - plane ): 



Im co 


co 


integration contour 


for large r 


saddle 

point 

4, 


T 


Re co 


FIGURE 37. Singularities in tu-plane for the BFKL equation with running QCD coupling. 


1. To draw contour through the saddle point ; 

2. To close contour on the zeros of the denominator . 


It turns out that the position of the saddle point u>sp = WL(r 0 ) and at large 
values of q 2 ( r ) it moves to the left. Therefore, for large r the rightmost singularity 
is the u; 0 (r 0 ) - the zero of the denominator. Closing contour on pole at u = co 0 we 
obtain Regge asymptotic cx oc e u ° y and all r- dependence enters only in the value 
of residue. 

For not very large r the saddle point contribution turns out to be dominant and 
leads to the Green function: 


G r {y ~yo,r,r 0 ) oc 


WL(y-yo) 

e 


( r — r n ) 2 

4 D ( y - y 0 ) 


(y-yo) 3 


(181) 


where aand D the running coupling constant ccs(r 0 ) should be taken at different 
scale: r ^ a instead of r 0 [49]. One can see that Eq. (181) gives a Regge-BFKL 
asymptotic of Eq. (19) but only at restricted values of y — yo < For larger 

“I 

y — y o the last term in exponent appears which has clear non - Regge behaviour 

[49]. 


C Summary. 

The BFKL Pomeron in pQCD gives an example of a different asymptotic be¬ 
haviour of the scattering amplitude than the Regge pole (“soft” Pomeron ). It is 
very instructive to understand that the space - time picture of the BFKL Pomeron 
is very similar to the “soft” Pomeron and a difference is only in the particular 
properties of the QCD cascade. 









For the running QCD coupling the BFKL asymptotics depend on the non per¬ 
turbative QCD initial condition which, in general, generates a “soft” Pomeron 
asymptotic behaviour at large values of virtuality. 

Certainly, the BFKL Pomeron is our window to non-perturbative regime and we 
need more experience especially with the BFKL asymptotic in the next-to-leading 
order approximation ( see Refs. [50] ) to make a reasonable educated guess about 
high energy asymptotic in non-perturbative QCD. 


X MY SEVERAL LAST WORDS. 


In this talk, I tried to explain what we have in mind when we are saying Pomeron. 
This is only an introduction in which I had to avoid many important details. I feel 
pity that I did not discuss the open theoretical questions and how the experimental 
measurements could help us to solve them. However, this is an interesting subject 
but certainly it cannot be considered in such an introductory talk. I also apologize 
that in the last section about the BFKL Pomeron, I was using mathematics too 
much, but I did this only to show you the only one example how a “soft” Pomeron 
can appear in QCD. 

I am very grateful to A. Santoro who invited me to his conference in Rio de 
Janeiro and without whom this talk will never be written. I thank all my brasilian 
friends:Francisco Caruso, Helio del Motta, Maria - Elena Pol and Ronald Shellard, 
whose hospitality recall me once more that I left a part of my heart in Rio. My 
special thanks go to E. Gotsman and U. Maor, who started with me re-analyze the 
situation in “soft” interaction more than five years ago and who spoke through my 
talk too, but with my full responsibility for all inconsistency here. 

I would like to acknowledge the kind hospitality of the Theory Group at DESY 
where this paper was completed. 
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